Depinning in a two-layer model of plastic flow 
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We study a model of two layers, each consisting of a d-dimensional elastic object driven over a random 
substrate, and mutually interacting through a viscous coupling. For this model, the mean-field theory (i.e. a 
fully connected model) predicts a transition from elastic depinning to hysteretic plastic depinning as disorder 
or viscous coupling is increased. A functional RG analysis shows that any small inter-layer viscous coupling 
destablizes the standard (decoupled) elastic depinning FRG fixed point for d < 4, while for d > 4 most aspects 
of the mean-field theory are recovered. A one-loop study at non-zero velocity indicates, for d < 4, coexistence 
of a moving state and a pinned state below the elastic depinning threshold, with hysteretic plastic depinning for 
periodic and non-periodic driven layers. A 2-loop analysis of quasi-statics unveils the possibility of more subtle 
effects, including a new universality class for non-periodic objects. We also study the model in d = 0, i.e. two 
coupled particles, and show that hysteresis does not always exist as the periodic steady state with coupled layers 
can be dynamically unstable. It is also proved that stable pinned configurations remain dynamically stable in 
presence of a viscous coupling in any dimension d. Moreover, the layer model for periodic objects is stable 
to an infinitesimal commensurate density coupling. Our work shows that a careful study of attractors in phase 
space and their basin of attraction is necessary to obtain a firm conclusion for dimensions d — 1, 2, 3. 



I. INTRODUCTION 

A. Overview 

Nonequilibrium transitions from stuck to moving states un- 
derlie the physics of a wide range of phenomena |T(j, from 
fracture and earthquake rupture J2IH0] to A ux fl° w m type-II 
superconductors pjlH] and sliding of charge density waves in 
metals 0, H, 0, llQfl - The rich collective nonequilibrium dy- 
namics of this broad range of phenomena can be modeled as 
an extended medium driven over quenched disorder. One can 
distinguish two main classes depending on whether the de- 
scription allows or not for plastic deformations of the medium. 
Within each class one may restrict to microscopic overdamped 
dynamics or allow for more complicated, e.g. inertial, effects. 

The first class of models, overdamped elastic media pulled 
by an applied force /, has been studied extensively. By defini- 
tion, the driven medium can be deformed by disorder but is not 
allowed to tear, and topological defects are excluded, the only 
degrees of freedom being elastic deformations. The question 
of applicability of this model to realistic situations is still de- 
bated in the static case (/ = 0), and even more so in the driven 
dynamics. The general expectation is that such a model is rel- 
evant to describe real systems in some range of lengthscales, a 
range which becomes broad (and potentially infinite, depend- 
ing on space dimension) in weak-disorder, strong-elasticity 
situations. Indeed one may conceive that, even if topolog- 
ical defects can be generated by the competition of elastic- 
ity, disorder and drive, they may remain bounded, and con- 
fined to shorter scales and thus unimportant for the effective 
large-scale description. This is known to happen in the stat- 
ics, e.g. for interfaces in random ferromagnets. Even when 
topological defects are relevant at large scale, the elastic de- 
scription may still apply at shorter scales. Hence the over- 
damped elastic model is a necessary first step to understand 
the collective dynamics in more complex situations. Many re- 
sults were obtained for this model, although some questions 



remain open: At zero temperature (T = 0) the elastic model 
exhibits a nonequilibrium phase transition from a pinned to 
a sliding state at a critical value f c of the driving force, first 
studied in the context of char ge d ensity waves 

El El- Start- 
ing from mean-field theory 11311 . an analogy with standard 
critical phenomena was developed, with the medium's mean 
velocity v acting as the order parameter, and a diverging corre- 
lation length iflL Il3fl . A functional extension of usual RG was 
developed to treat quenched disorder and obtain the rough- 
ness and dynamical exponents at the threshold v = + to 
1 -loop accuracy ITil [l5ll . Extensions at non-zero v empha- 
sized the differences with standard critical phenomena llql . It 
was shown that a two-loop Functional RG (FRG) approach is 
necessary to fully describe the difference between statics and 
v = + quasi-static depinning lfl7i [l8i [l9ll and to reach sat- 
isfactory agreement with numerical simulations l2(il l2ll l22ll . 
Universality classes were identified, which are distinguished, 
for example, by the range of interactions or by the periodicity 
(or nonperiodicity) of the pinning forces. A key feature of the 
overdamped elastic model is that for one component displace- 
ments N = 1, i.e. interfaces, the sliding state is unique, the 
v(f) curve is single-valued, and no hysteresis can occur in the 
moving state at v > 0. This property, based on Middleton's 
theorem l23ll . which also leads to simplifications lfl8ll in the 
FRG description for N = 1 is not expected to hold for N > 1. 
As a result, the understanding of the > 1 depinning tran- 
sition for e.g. lines or vortex lattices, is still not satisfactory 
despite some attempts I24ll25il26ll . Furthermore, there is a sec- 
ond type of universality classes for depinning (e.g. anisotropic 
depinning) which does not obey the so-called statistical tilt 
symmetry (or rotation symmetry) and where non-linear terms 
become relevant (e.g. Kardar-Parisi-Zhang (KPZ) like terms) 
ll27ll . Despite efforts lf28[ |29ll . a complete theory for this class 
is still lacking, and even the value of the upper critical dimen- 
sion is a matter of debate. The question of non-linear terms 
may be of importance to experiments of contact line depinning 
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and cracks ||3(X|3U 13211 . Away from depinning, well into the 
uniformly sliding state at v > 0, it was found that the dynam- 



ics can be surprisingly rich 
ponent periodic objects |[34] 



13311 . especially for N > 1 com- 
35h . New terms can be generated 
in the equation of motion, a linear convective term, a static 
random-force term, and a host of pos sible non-linear, KPZ- 

type terms OS USUI HIS SHE EI B- For N > 1, 
a distinct, "moving glass" fixed point was found in the FRG, 
with persistence of transverse order and transverse pinning, 
leading to the prediction of a moving Bragg-glass and a mov- 
ing smectic state l38ll42l l43ll . In both states the flow is orga- 
nized in static-like channels, in a layered fashion. Extensions 
to correlated disorder was studied, and a moving Bose-glass 
state predicted l45l l46ll . Although clear evidence of these ef- 
fects were found in numerics and experiments ll44tl . no sys- 
tematic study of finite-size corrections was made. Since the 
simultaneous analytical treatment of all the terms allowed by 
symmetry within a FRG approach is a problem of formidable 
complexity |43j, even a fully consistent theory of the elas- 
tic flow at large velocity is still lacking. Hence the question 
of which moving state is stable in the thermodynamic limit 
is still open. Finally, once the elastic system is understood, 
one may hope to construct arguments for or against stability 
of the elastic flow to defects. These however are even more 
delicate than in the statics, where the stability of the Bragg 
glass was debated, and the validity of the driven elastic model 
has only been assessed qualitatively l42ll43T l47l l48ll . Hence, 
as one can see, despite being well studied, the overdamped 
elastic model is still far from being understood. Extensions to 
include inertial effects and stress overshoots 1 have also been 
considered |49ll , but much work remains to be done. 

There are many experimental situations where the elastic 
medium model seems insufficient and one needs to take into 
account plastic deformations, as e.g. topological defects in 
the medium. In a wide class of experiments strong disorder 
yields large deformations of the driven medium that make 
a strictly elastic model of the extended structure inapplica- 
ble mumm. 

In contrast, the medium tears as topo- 
logical defects are constantly generated and healed by the in- 
terplay of drive, disorder and interactions 11531 154 p5l Da 1571 
l58l l59l l60ll . At slow average flow rates the dynamics near 
depinning is spatially and temporally inhomogeneous, with 
coexistence of pinned and sliding degrees of freedom. The 
depinning transition may become discontinuous (first order), 
possibly with a macroscopic hysteresis and switching between 
pinned and sliding states J6ll l62tl . Experiments on charge 
density waves show that varying the temperature leads to a 
transition from continuous depinning to hysteretic depinning 
with sharp switching between pinned and sliding states l6lll . 
Whether such phase slip effects occur in the bulk or only at 
the contacts ]63[], remains to be clarified. Related slip effects 
or plastic behavior have been proposed to explain the com- 



1 Note that in the mean-field limit the stress-overshoot model and the model 
studied here (see below) are identical, aside from the fact that in the crack 
model nonperiodic disorder has been considered. 



plex dynamics of many other dissipative systems, including 
vortex arrays in type-II superconductors. Lorentz microscopy 
images of driven vortex arrays in irradiated thin films of Nio- 
bium show vortex rivers flowing past each other at the bound- 
aries of pinned regions of the lattice 115 311 . Scanning tunneling 
microscopy, which can resolve individual vortices at high den- 
sity, also reveals the evolution of the vortex dynamics with 
disorder strength l64ll . There too, there are edge contami- 
nation effects, and they may be responsible for the coexis- 
tence of a metastable disordered phase and a stable ordered 
phase ll65ll66ll . It is clear that more work is needed to under- 
stand the rich dynamics of driven systems in experiments. 

It was ubiquitously found in numerical studies of interact- 
ing particles driven on a random substrate at T = (away 
from the weak-disorder limit) that near the onset of mean slid- 
ing the motion occurs along filamentary channels or rivers that 
are determined by the spatial disorder of the random medium. 
Such channels are preferentially aligned along the direction 
of mean motion, but can exhibit large transverse excursions. 
At higher mean-flow rates the rivers coalesce into a more 
coherent structure that eventually results in a uniform flow. 
Hence the plastic flow takes, at a qualitative level, a vari- 
ety of forms with increasing correlations: (i) filamentary flow 
with a single well-defined channel or several uncoupled chan- 
nels 115 2[ l67l l68l l69l 17011 to coupled or synchronized chan- 
nels, to a layered smectic type structure to a moving lattice 
which may or not still contain frozen or moving dislocations 
|H HI UlEU HUH. While one may hope that at large 
velocity, where the effective disorder is smaller, the flow is 
closer to the one described by an elastic model, it is clear that 
one needs to take into account plastic deformations to describe 
these various regimes. 

The theoretical understanding of the dynamics of such plas- 
tic systems is much less developed than that of driven elastic 
media. It is not even clear how to characterize the various 
moving states which are observed by some order parameter, 
and to properly define steady states and their large-size limit. 
One can measure the distribution of time-averaged velocities 
P(v) of the individual particles. A non-trivial P(v) exists for 
instance in the filamentary regime where some particle seem 
permamently pinned while others are moving along channels. 
In small systems with periodic boundary conditions a peri- 
odic steady state is observed near the threshold with a non- 
trivial P(v). Whether this feature persists in the infinite sys- 
tem limit, and how it depends e.g. on the geometry and aspect 
ratio of the sample, is not known. As was recently pointed 
out l67ll . it is fruitful to apply tools and ideas from the the- 
ory of dynamical systems and chaos. It was found that upon 
increasing /, the system undergoes a transition from periodic 
to a fully chaotic flow with positive Lyapunov exponents and 
a non-trivial attractor. The dimension of this attractor, which 
is low, may also provide a tool to characterize the phases of 
plastic flow. These ideas remain to be explored, in particular 
whether the elastic flow could exhibit some chaotic regime. 
At larger drive, P(v) becomes more peaked around a single 
velocity and some degree of spatial coherence in the phase 
accross the layers may arises. Whether P(v) eventually be- 
comes a delta function in the large-size limit, and whether 
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the phase-coherence lengths diverge or not, has not been sys- 
temetically studied numerically. There has been some efforts 
to use numerical simulations to correlate the spatial and tem- 
poral structure of the dynamics with the shape of the macro- 
scopic response l60h . for instance the IV characteristics due 
to flux flow in a type-II superconductors. A number of mean- 
field models of driven extended systems with locally under- 
damped relaxation or phase slip s have been proposed in the 

literature mmizlHlzlUflSlziHiilsS]- Whether or 
not these dynamical models exhibit truly collective behavior 
and universality in finite dimensions remains an open ques- 
tion which motivated this work, as discussed below. A model 
which attempts to describe filamentary flow away from mean 
field was proposed in l82ll . 

B. Layered Model 

Given the difficulty in describing topological defects, a sim- 
pler approach consists in considering layers such that defor- 
mations within a layer are only elastic. Since the relative dis- 
placements between layers can be arbitrarily large, inter-layer 
plastic deformations are allowed. Whether they occur or not 
depend on the interaction between the layers. This approach 
was successful to treat disorder in the statics, where it lead to 
solvable limits for e.g. the Bragg glass phase lf83"ll84ll . the de- 
coupling transition for magnetically coupled superconductors 
JUt]. It is also studied to describe interacting quantum systems 
such as the sliding Luttinger liquid ll86ll . Recently a similar 
strategy was applied to describe plastic flow and depinning 
(see Refs. l87l l88tl for a review), and coupling phenomena 
in the driven dynamics ll76ll . There it is even more natural 
since the flow naturally tends to be along layers (which can be 
channels) in the direction of the applied force. In one version 
of the model, introduced by one of us and collaborators, the 
layers are only viscously coupled in at least one of the direc- 
tions transverse to the mean motion. Although this is a much 



simplified version of plastic flow, for instance there are no 
convective terms in the equation of motion, it is motivated by 
the moving smectic phase in driven vortex lattices mentioned 
above. It incorporates elastic responses to compressional de- 
formations and allows for local slips of neighboring degrees of 
freedom due to shear deformations. It is also relevant to exper- 
iments on driven superconducting vortices in narrow channels 
and other controlled geometries l89l \9(X \9 ll l92l 19311 . One pos- 
sible realization is a a layered (e.g. high T c ) superconductor 
when the vortices are aligned with the magnetic field within 
the CU2O, a6-plane layers and move along these layers under 
a c-axis current. In the limit where the intrinsic pinning po- 
tential from the CU2O planes is strong compared to the weak 
isotropic disorder from point impurities, the vortex dynam- 
ics may be modeled in terms of 2d elastic layers or "chan- 
nels" coupled viscously along the c axis. The fact that only 
the viscous coupling between layers is retained makes it more 
tractable. It is expected to be valid in situations where the 
commensurate density-density interlayer interaction, studied 
in Ref. ||7(| , which couples the displacements in each layer, 
can be neglected. The general case can be defined as fol- 
lows: Consider a d = dii + d± -dimensional medium com- 
posed of elastic dii -dimensional channels coupled via viscous 
interactions in the remaining d± directions. The medium is 
driven by a uniform force f applied along one of the direc- 
tions in the dii -dimensional channels. Here we only consider 
the dynamics of a scalar displacement field u(ri|,r_i_,i) de- 
scribing deformations in the direction of the driving force at 
position r = (rn , rj_), with r, rn and rj_ vectors in d, dy 
and d± dimensions, respectively. To index the channels one 
discretizes spatial coordinates in the direction normal to the 
layers (rj_ — > r n , where r n denotes the n-th layer) and let 
r 11 = x. The dynamics of the displacement u n (x, t) of each 
degree of freedom is governed by the equation, 
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ryd t u n (x) — / K(x — x'){u n {x) — u n (x')) + ^ (x) - ii n (x)] + f + F n (u n (x), x) 



(1.1) 



This is the M-layer model. Among the various universality 
classes of disorder, the one of most interest here is the random 
periodic class where the pinning force has the form: 

F n (u n ,x) = hiY{u n {x) - /%) , (1.2) 

with Y(u) a periodic function. The pinning strengths h % n 
are independent random variables distributed with probability 
p{h z n ) and (3 l n are random phases uniformly and independently 
distributed in [0, 1). This models the dynamics of driven pe- 
riodic media, such as vortex lattices, charge density waves, or 
Wigner crystals. In these systems substrate impurities couple 
to the density of the lattice which, in the absence of in-layer 
topological defects, has the periodicity of the ordered lattice. 



As a result, the pinning force contains periodic components 
at all reciprocal lattice vectors 15l|94j]. In the bare model one 
can retain the lowest Fourier components only, since as is well 
known from FRG studies of statics and depinning, all Fourier 
components are generated by coarse-graining and should be 
included to describe the properly renormalized disorder cor- 
relator. Such a correlator develops cusp-like singularities at 
large scales that control the dynamics. The other type of disor- 
der, the non-periodic or random-manifold class, which at the 
elastic depinning was shown to give rise to a single universal- 
ity class encompassing both random bond (i.e. short-ranged) 
and random field (i.e. long-ranged) disorder, will also be stud- 
ied. This is done by choosing a non-periodic correlator for 
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the random forces F n (u, x). Physical realizations are less ob- 
vious since the above velocity coupling is local in x space 
only while a realistic coupling e.g. between two directed poly- 
mers would also depend on the field u. Two possible realiza- 
tions are: (i) manifolds with internal disorder, as studied in 
Ref. ll95i 19611 . (ii) periodic systems for which the correlation 
length of the disorder r/ is small compared to the lattice spac- 
ing a. Then the two scales for pinning, the Larkin length R c , 
and R a for the decay of translational order, can be very dif- 
ferent, and it is known that for scales R c <C L <C R a all har- 
monics of the disorder correlator are important and the system 
behaves effectively as a random manifold within this range of 
scales JH HH H3]. Hence, below we also consider the non- 
periodic or random-manifold model and discuss the different 
behaviors in the two cases. 

C. Aim and outline of the paper: Two-layer model 

The layered model ( II. lb with viscous couplings was proposed 
as a generic coarse-grained model representative of a class of 
dissipative driven systems with strong disorder that encom- 
passes many of the models considered in the literature. It 
was solved within mean-field theory and shown to predict a 
qualitative change from continuous to discontinuous and hys- 
teretic dynamics as a function of disorder strength, consistent 
with experimental observations in a variety of systems l8(ill . 
It has also been studied numerically in finite dimensions. The 
numerical studies show evidence of hysteresis in 2+1 dimen- 
sions above a critical value of the interlay er coupling. Hys- 
teresis was not clearly evident, however, in 1+1 dimensions 
nor for the two-layer model studied below, although it could 
also not be conclusively ruled out on the basis of finite-size 
scaling rfsTTl - 

The aim of this paper is to go beyong the mean-field treat- 
ment of model ( II. j} and explore using functional RG whether 
hysteretic dynamics also occurs and whether universal fea- 
tures emerge in low dimensions where one usually does not 
expect the mean-field approximation to be accurate. Since 
generalization to M layers is straightforward and not expected 
to bring important qualitative changes, we study in detail 
the technically simpler case of two viscously coupled layers 
M = 2. We start by recalling in Section [TT] the main fea- 
tures of the mean-field solution so as to provide a basis for 
comparison. In Section [HI] we study the M = 2 model first 
by direct perturbation theory and next using 1-loop FRG. We 
prove that the elastic single-layer (i.e. decoupled layer) quasi- 
static v = + depinning fixed point is always unstable to a 
small viscous inter-layer coupling. A partial one-loop analy- 
sis at v > shows the generic co-existence of a pinned and a 
moving state below the single-layer depinning threshold. The 
resulting v(f) curves show similarities with the mean-field 
ones, and in some regimes the agreement can even be made 
quantitative. We estimate the velocity v c at which the v(f) 
curve becomes vertical (and a jump may occur in the fixed 
force ensemble). A key feature of the one-loop study is that 
the inter-layer viscous coupling 712 is not corrected by disor- 
der. To determine whether this is maintained to higher order, 
we carry in Section[lV]the FRG to two loops, near the uncou- 



pled elastic quasi-static depinning. The analysis gives a strong 
correction as compared to one loop in the non-periodic class, 
i.e. random manifolds, with a new universality class for plas- 
tic depinning whose z exponent is computed in Section lTV Bl 
The analysis in the periodic case is quite subtle and presented 
in Section lTV CI Finally, to get a better understanding of pos- 
sible behaviors and of the connection between dynamical hys- 
teresis and attractors in phase space in simpler cases, we study 
in Section [V] two d = toy models of two viscously cou- 
pled particles on, respectively, a smooth (Section [V Aj and a 
discontinuous (Section TV B| ) force landscape. Finally, in Sec- 
tion lV Cl the main results are summarized, and extensions and 
future directions discussed. In particular it is proved that a 
small interlayer commensurate coupling is irrelevant at depin- 
ning. Since such a term is always present in real systems, this 
shows that the viscous model is consistent. The Appendices 
contain the details of the two-loop calculation and a proof that 
stable static configurations where decoupled layers are inde- 
pendently pinned remain dynamically stable in presence of the 
inter-layer viscous coupling. 

Let us now define the two-layer model studied here, and 
fix notations. We consider the overdamped dynamics of two 
layers coupled by a viscous coupling in a random poten- 
tial. Each layer is an elastic system parameterized by a one- 
component (N = 1) displacement field u l (x), also denoted 
u x , for i = 1, 2, or u x t to indicate explicitly the dependence 
on time. The equation of motion of one layer is 

Tb< t = m{ul tt -u 1 Xtt )+cV 2 u 1 Xtt + F 1 (u 1 Xtt ) + f , (1.3) 

where 70 is the in-layer friction coefficient. Hence, in addition 
to elastic intra-layer restoring forces (elastic coefficient c) and 
the quenched random pinning force, one layer is also pulled 
by the other layer through a velocity (or viscous) coupling i]q. 
Here we focus on the case of uncorrected disorder in each 
layer, and denote the second cumulant of the pinning forces 
by 

F i (x,u)Fj{x , ,u')=6 ij S d (x-x')A (u-u') . (1.4) 

The equation of motion for the system of two layers driven by 
an external force / can then be written as: 

( in 712 \ d fulA_ 2 ful t \ fF 1 (x 7 u 1 xt ) + f\ 
{112 122) dt [ulj - CV {ulJ + {F^x,ul t ) + f) 

(1.5) 

The bare values for the friction matrix are 

7n = 722 = 7o + Vo (1-6) 
712 = -Vo ■ (1-7) 

II. MEAN-FIELD THEORY 

To set up the mean-field theory for the multi-layered model, 
it is convenient to discretize space in both the transverse and 
longitudinal directions, using integer vectors t, m for the dn- 
dimensional intra-layer index. The local displacement along 
the direction of motion at time t is u\(t), with i = 1, . . . , M 
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the layer index and I — 1, . . . , N labeling the degrees of free- 
dom within each layer. Its dynamics is governed by the equa- 
tion (in this section we drop the subscript '0' on the bare fric- 
tions), 

7tiS(t) = £ K lm [u\ n {t) - 4(t)] + ]T thMW - 



(2.1) 



where Y(u) = Y(u + n) is a periodic function and K( m and 
%j have constant values c and 77, respectively, for nearest- 
neighbor pairs and vanish otherwise. The random pinning 
strengths h\ are chosen independently with probability dis- 
tribution p(H\) and the random phases (3\ are distributed uni- 
formly and independently in [0,1). 

A. Fully connected mean-field theory 

One mean-field approximation is obtained by assuming that 
all sites are coupled with uniform strength, both within each 
layer and across the layers, i.e., K( m = c/N for all I and m 
and rjij = rj/M for all i and j. The mean displacement and 
velocity are given by 



u(t) 



— YY4(t) 

i i 



NM 



(2.2) 
(2.3) 



and we look for solutions moving with a uniform velocity so 
that (up to a choice of the origin of time) 



l(t) = vt 



(2.4) 



Since the displacements are coupled only through the mean 
fields, they can be indexed by their disorder parameters (3 
and h, rather than by the spatial indices i, i, i.e., u\(t) — > 
u(t; (3, h). The mean-field dynamics is governed by the equa- 
tion 

(7 + r])u(t; (3, h) = c(vt -u)+f + r]V + hY(u - (3) (2.5) 

that must be solved with the self-consistency condition that 
determines the mean field, 



{u(t;/3, h) - vt)f3. h = , 



(2.6) 



where (■■■)/3.h = f d(3 J dh...p(h) denotes the average over 
disorder. 

The long-time steady-state solution to Eq. ( 12. 5i can be writ- 
ten as 



u(t; f3,h) = vt + u , 



with 



(7 + = —cu + / — jv + hY{u + vt — (3) 



(2.7) 



(2.8) 



to be solved with the condition (ii)p.h = 0. It is apparent from 
Eq. ( 12.81 ) that u is a periodic function of time with period 1/v 



and depends on time and phase (3 only through the combina- 
tion u = u(t — v//3; h). This will allow us to carry out the 
average over (3 by averaging over time. 

We display below the solution for a parabolic scalloped pin- 
ning potential, corresponding to a piecewise linear pinning 
force with jumps of size h at the boundaries of each period, 



Y(u) 



n < u < n + 1 



(2.9) 



with n an integer. The mean-field equation ( 12.51 ) is for- 
mally identical to the mean-field equation for a purely elastic 
medium, with friction coefficient T = 7 + rj and an effective 
drive F = f + i]v. The solution of the mean-field equation 
for a scalloped pinning potential and r\ = was obtained by 
Narayan and Fisher lfl4f1 and is easily adapted to our case. The 
solution for finite 77 is described in Ref. 1801 18811 and will be 
summarized here for completeness. 

The pinning force has a jump discontinuity h at the end of 
each period. The displacement u is continuous across neigh- 
boring periods, but the local velocity u has jumps of size h/T 
at tj + n/v. The solution of Eq. ( |2.8t for tj((3) + n/v < t < 
tj(J3) + (n + l)/v is 



Ae 



_ xt h(/3-vt) + F + h(n+ 1/2) cv 



(2.10) 



TX TX 2 ' 

where A = (c + h)/T and tj((3) is the "jump time". The 
constant A(f3) and the jump time tj((3) are determined by 
requiring 

u(t = tj + n/v) = n + (3 , (2.11) 
u(t = tj + {n + \)/v) = n + 1 + ft . (2.12) 

It is important to appreciate a crucial difference between the 
mean-field theory of the purely viscous model (c = 0, or 
dii = 0) discussed in Refs. J77|,|78|] and the mean-field theory 
of the model considered here that includes additional elastic 
couplings within the channels. In the purely viscous case, 
each degree of freedom is coupled only to the local veloci- 
ties (which exert an additional effective driving force) and can 
slide with its own period. In contrast, when c ^ each de- 
gree of freedom couples to the average displacements via a 
spring-type interaction that forces all periods to be the same, 
independent of h. 

After inserting A and tj obtained from the solution of 
Eqs. ( 12.1 It and d2.12| i in Eq. (12. 1 0b , it is straightforward to 
impose the self-consistency condition as 

u(t-0/v;h)) =v( [ " dtu(t~/3/v;h)) =0. 

/ 0,h \ Jt, + !± ' h 

(2.13) 

This yields an implicit solution for the mean velocity as 

F(v) = f(v) +r,v = f c +Tv[l- Af(c)] 

h 2 1 \ 

(2.14) 



.(c + h) e x / v - l/h 
with f c the threshold force for the onset of uniform sliding, 

h 2 



fc (2(0 +h))h 



(2.15) 
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FIG. 2.1: The mean-field velocity «(/) as a function of / near thresh- 
old, as given by Eq. J2.17b . The slope of the linear function diverges 
at r\ — rj c (c) and is negative for 77 > T) c (c). 



and 



M(c) 



(c+/i) 5 



(2.16) 



The threshold force depends only on the elastic coupling c, 
but not on the viscous coupling 77. This follows because in 
mean field the viscous coupling becomes effective only when 
the system is moving as a whole, while away from mean field 
one expects additional fluctuation effects. 

Near threshold the term on the second line of Eq. ( 12.14b 
gives contributions of order ~ e -1 / 1 ' and can be neglected. It 
is then easy to invert Eq. ( 12.141 ) to obtain v(f), with the result 



(/ - .0 



7 — (7 + r))M(c) 



(2.17) 



The mean velocity vanishes linearly for / — ► /+ , with a 
MF exponent /3mf = 1 which is generic for discontinuous 
pinning forces fTill . When 77 = the slope of the linear curve 
is always positive as M(c) < 1. The slope diverges, however, 
at a critical value of 77, 



Vc(c) = 7 



1 



M(c) 



1 



(2.18) 



and becomes negative for 77 > ?7 C , as shown in Fig. 12. II For 
T] > rjc(c) the velocity curve is multivalued, yielding hys- 
teretic behavior. 

The phase diagram and typical velocity-force curves are 
shown in Figs. I2.2l and 12 . 3 1 for p(h) = S(h — ho). The finite 
long-time elasticity (c 7^ 0) guarantees that the behavior is in- 
dependent of the shape of the pinning-force distribution p(h). 
The phase diagram for p(h) = e~ h was shown in i80ll and has 
the same form as the one shown here. The point (r/ c , f c ) is 
a tricritical point separating single-valued from multi-valued 
velocity curves. For 77 < 7y c , a continuous depinning transition 
at f c separates a pinned state from a sliding state with unique 
velocity. A question addressed below is whether i] c remains 
non-zero in finite dimension and if so, whether the depinning 
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FIG. 2.2: Phase diagram obtained from the fully-connected mean- 
field solution for p(h) = S(h — 1), for 7 = 1 and c = 1, corre- 
sponding to 77 c = 3. There is a critical point at (rj c , f c ) separating 
continuous from discontinuous depinning. 
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FIG. 2.3: Typical velocity-force curves obtained from the fully- 
connected mean-field solution for p(h) — 8(h — 1), c = 1 and 
rjc — 3. For 77 = (dashed curve) and 77 = 1 (solid line) the system 
depins continuously at f c . At 77 = ?7 C = 3 (dotted line) the slope di- 
verges at threshold. For 77 = 12 (solid line) the velocity-force curve 
is multivalued. This corresponds to a hysteretic depinning transition 
as the system depins at f£ when the force is ramped up from zero 
and repins at the lower value when the force is ramped down 
from the sliding state. 



transition for r/ < r} c is in the same universality class as the 
depinning of an elastic medium (77 = 0) lfl3ll . 

In our mean-field example, the linear response diverges at 
77c as ^(77 = f] c ) ~ l/ln(/ — f c ). For 77 > ?7 C the solution 
is multivalued. In this case when the force is ramped up from 
zero the system depins at = f c . When the force is de- 
creased from a value above f c the system gets stuck at the 
lower value ff, yielding hysteretic v(f) curves. 

For 77 > i] c the mean velocity has a jump discontinuity. The 



7 



value v c of this jump is given by the solution of 

V OV ) V=V C 



v(y+r\ ) 



(2.19) 



where f(v) is given by Eq. (12.14b . An explicit solution for 
the jump v c can be obtained for the case of a sharp disorder 
distribution p(h) = 5(h—ho). In this case the condition (12. 19t 
for the jump becomes 



V - Vc 



(X/v) 



7 + V 4sinli (\/{2v)) 
For rj ^> 7, ?7 C this gives 



/in 



'/ 



2\/377 V 7 + 57- 



(2.20) 



(2.21) 



Finally, we note that the mean-field theory for a smooth peri- 
odic pinning potential gives qualitatively the same phase dia- 
gram, although with mean-field exponent 0mf = 1/2. 

The fully connected mean-field theory discussed here for 
the layered visco-elastic model is formally identical to the 
mean-field limit of a model of crack propagation with stress 
overshoot studied by Schwarz and Fisher lf7ll p72l l98tl . al- 
though the crack model contains random force disorder in- 
stead of the periodic disorder considered here l87ll88ll . 

B. Self-consistent single-layer approximation 

An alternative, "partial" mean-field approximation treats only 
one direction of space using mean field, and reduces the prob- 
lem to an effective single-layer model. It is obtained by 
assuming uniform, i.e. infinite-range couplings of strength 
Vij = Tl/M across the layers for each in-layer site I. The 
corresponding mean field is given by 



(2.22) 



In the thermodynamic limit of an infinite number of layers, 
assuming the system is self-averaging, the mean field ve will 
not depend on I and this label can be dropped. The mean-field 
dynamics is then described by the equation 

( 7 + v)4(t) = E - «<(*)] 

m 

+f + T)v + hiY(4(t)-f3i), (2.23) 

which must be solved with the condition {vfyp^h = v - It is 
illuminating to rewrite Eq. ( 12.231 ) by replacing the discrete in- 
layer index I by the original continuum variable x, 

( 7 + V )u* = cV 2 < t + / + V v + KY(ul t -PI), (2.24) 



l x,t 



to be solved with the self-consistency condition -j ^ 
v. It is apparent that Eq. ( 12.241 ) describes the dynamics of 
M identical elastic layers coupled only through the mean 
field v. Each layer is a dissipative elastic medium of friction 
r = 7 + 17, driven by a force F = f + rjv. For 77 = 



• (f-fj 



fc 



f+ J] V 



FIG. 2.4: The mean velocity (7 + rf)v plotted as a function of F = 
f+r]V. When f+rjv is used as the independent variable, the velocity- 
force characteristic has the same functional form as that of a single 
elastic layer that depins at a threshold f c with v ~ (/ — /c)' 3 and 
/3<las/^/+. 



the layers are decoupled, with r = 7 and F = f. The 
velocity-force curve v s i(f) of one deco uple d layer has been 
studied in details H [3 H H |H Each layer is 

pinned with v = for / < f c . It depins at / = f c and 
slides for f > f c with mean velocity v s \(f) = Q{f)/"t, 
and G(f) ~ (/ - f c f as / /+ and (3 < 1 a'criti- 
cal exponent that depends only on the system's dimension- 
ality. It is clear from the form of Eq. ( 12.231 ) that the velocity- 
force characteristic of the coupled layers has the same func- 
tional form as that of an individual layer, with the replacement 
/ -> F, i.e, v(F) = Q{F)/T. A sketch of this velocity- 
force characteristic is shown in Fig. 12.41 The velocity-force 
characteristic v(f) of the coupled layers can then be obtained 
simply by performing a shift in the independent variable in 
the known result for a single layer. The result is shown in 
Fig. ES Near threshold v(f) ~ (/ + rjV - fc) 13 , with 
P= l-(4-d||)/6 + C[(4-d||) 2 ] < 1 and v(f) will be mul- 
tivalued for every finite value of rj, yielding a hysteretic depin- 
ning transition 118011 . The hysteresis for any 77 > 0, for d < 4, 
obtained here in the approximation of a global transverse cou- 
pling, will be confirmed below within a one-loop FRG analy- 
sis which incorporates inter-layer fluctuations neglected here. 
In both cases it is a consequence of the non-trivial renormal- 
ization of 7 within a single layer, responsible for j3 < 1 and 
z > 2 for elastic depinning. 



As pointed out in Ref. 18011 . the self-consistent single-layer 
approximation, with uniform couplings across the layers, is 
equivalent to a model of charge density waves (CDWs) that 
incorporates the coupling of the CDW to normal carriers via 
the addition of a global velocity c oupling to the equation of 
motion for the phase ifTl [lQll [Tol . 

Finally, another "partial" mean-field theory is obtained by 
assuming uniform couplings of strength K( m = K/N for ev- 
ery £, m within each layer. This model will be discussed else- 
where. The two-particle toy model described in Section V 
corresponds to the K = limit of this mean-field theory. 
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FIG. 2.5: The velocity-force curve for finite r\ can be obtained fro the 
single-layer curve of Fig. l2.4l corresponding to r\ — (dashed line) 
by a change of the independent variable. Since /3 < 1, the resulting 
v(f) will be multivalued for any finite rj. 



III. FUNCTIONAL RG TO ONE LOOP 

To go beyond mean field we now develop a Functional RG 
approach. 



A. Perturbation theory and length scales 

Consider model ( 11. 5\ driven by a force / and assume that it 
reaches a time-translational invariant steady state (e.g. with 
periodic boundary conditions for each layer). There are two 
modes: 



u+ = (u 1 +u 2 )/2 



(3.1) 



For a system of finite size L, because of fluctuations in the 
pinning force, the velocity in each layer will be different. 
However this effect should disappear in the infinite-L limit, 
and can be supressed using appropriate boundary conditions. 
Hence we define v to be the velocity of the center of mass 
V = u + , perform the shift to the comoving frame, 



vt 



(3.2) 



and immediately drop the hat. We can now write the dynam- 
ical action associated to the resulting equation of motion (i.e. 
CT3) shifted): 



S[u, u] 



U°2 7?! 



x,t,V 



At 

Mult 



l x,t' 





-cV J 

v(t - 



"x,t 
l x,t 



if 



7» 



t')) 



M<t 



+ v(t-t')) 



(3.3) 



The subscript indicates that these are quantities for the bare 
model (and it is often dropped in the following). The matrix 
of friction coefficients is diagonal in the basis (13. U and we de- 
note the frictions associated to the center of mass and relative 
motion as 



7+ = 7n + 7i2 



7- 



7n — 7i2 



(3.4) 



The bare values are 7+ = 70 and 7° — 70 + 2r/o. The bare 
response functions, i.e. those in the absence of disorder, read: 



n k,t 



'fet«_fc n 



R 



11 _ D 22 



k,t 



Kit = m 



-tk z h 



+ 



12 



r?2i 

n k.t 



Q(t) 



27+ 

e -*fc 2 /7 



+ 



20 



+ 



27- 
e -tfc 2 /7- 
27- 



(3.5) 

(3.6) 
(3.7) 



The case of a single layer is reproduced upon setting 770 = 0, 
or equivalently 7_ = 7+ (then Rj^ — > Rut, the standard 
single-layer response function, and R]^ t — > 0). 



Writing S = So + Si n t where Si n t contains only the dis- 
order, i.e. the second line in (13.3b . the effective action T[u] of 
the system can be computed perturbatively in the disorder: 

T[u,u] = S [u,u] + (S- mt [u,u})s + ^(S int [u, u] 2 ) c So H 

(3.8) 

In the average over only 1 -particle irreducible (1PI) di- 
agrams (i.e. containing loops) are kept. The quadratic part 
of the effective action yields the exact disorder-averaged re- 
sponse and correlation functions: 



,X -q,t-t' 



U qtU_ qt , 



S 2 T[u,u] 



5ui qt ,8u 3 qt 

~ ^q,t-t' 



(3.9) 



u—u—O 
k 



-qti uu qt 2 



u—u—O 



Both functions are symmetric in i, j and in q. The effective 
action has a complicated form but contains terms similar to 
those in the above action (13.3b with renormalized (i.e. "cor- 
rected") values for the friction matrix 7^ , and the second cu- 
mulant of disorder A(u). The elastic term is unrenormalized 
(i.e. the zero frequency part of the uu term in T is the same as 
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in So) thanks to the statistical tilt symmetry 2 which holds in- 
dependently in each layer. Other terms are generated in pertur- 
bation theory, such as higher disorder cumulants, higher fre- 
quency corrections to the self energy or non-linear terms such 
as J u(dtu) 2 . In each case their relevance should be assessed 
carefully. These terms are usually irrelevant near d = 4. A 
simplifying feature is that the coupling between the layers is 
purely dynamical. Therefore the static part of the theory (i.e. 
the O-frequency part of the effective action) consists of two 
decoupled static layers. This implies, among others, that no 
outer-diagonal elements of the disorder correlator are gener- 
ated in perturbation theory. 

Let us now examine perturbation theory and power count- 
ing. The effective action contains the term: 



-/ 



f = f -lov + 5f(v) (3.10) 



where Sf(v) contains all corrections due to disorder. On av- 
erage these are the same for each layer, and depend on v. 
The equation of motion is obtained from the condition / = 
equivalent to (u xt ) = (in shifted variables). From ( 13 .8b one 
finds that to lowest order in A (i.e. to one loop) the corrections 
to friction and force are: 



#7i 







*7n = - / / dTTA"(«r)i# 

Jq JO 



Sf 



dr A'(vt)RII , 



(3.11) 
(3.12) 

(3.13) 



qJO 



where the index is implicit if one studies perturbation theory 
on the bare action. The correction to the disorder A is of order 
A 2 and, at v = 0, is identical to the one for a single-layer 
model, while at finite v it has a complicated expression (even 
in the single-layer case, as given in lUql not displayed here. As 
is well known, for v = + , A(u) acquires a cusp for scales 
larger than the Larkin length L c . 

Before obtaining the 1-loop FRG equations let us make 
some general qualitative comments on the stability of the 1- 
layer elastic quasi-static depinning to the viscous inter-layer 
coupling. The absence of 1-loop corrections to 712 implies 
that to this order 712 = 7i 2 = —rjo- Consider quasi-static 
depinning v = + . Then one finds 



57i 



~7nA"(0 H 



<f 



(3.14) 



where a UV cutoff is implicit everywhere. This is the same 
correction as for the single-layer problem (i.e. for v = 
it does not depend on 712); hence under coarse graining 
711 is reduced compared to its bare value (above the Larkin 



2 The invariance of the non-linear (i.e. disorder) terms of 5 under u L — ► 
utf + <j> l (x) for two arbitrary static functions <j> l (x), i = 1,2, should 
persist for I\ 



scale L c , A"(0 + ) is strictly positive). The intra-layer fric- 
tion 711 < 7°j = 70 + 770 remains finite and non-zero 
for d > 4 (where the above integral converge at small q) 
while for d < 4 it becomes dependent on the system size 
L, 711 ~ (70 + Vo){L / L c ) z ~ 2 , z < 2 being the single-layer 
dynamical exponent for elastic depinning in d < 4. Since 712 
is uncorrected (it is negative) and 711 is reduced, it is clear 
that the friction coefficient of the center of mass of the system 
7 + = 711 + 712 may become negative at some scale, denoted 
L p \. When this occurs the fixed point of elastic quasi-static 
depinning becomes unstable (and inconsistent). This always 
occurs for d < 4, but only for 770 larger than a critical value rj c 
for d > 4. The qualitative picture is then as follows: 
(i) d > 4: To lowest order the equation of motion reads: 

[«(7o + m) - Vo] v = f-fc + 0(v\) (3.15) 
K = W7n = l + A"(0+) f^T, (3.16) 



where we denote by < k < 1 the usual reduction factor 
in friction in the single-layer problem. Elastic quasi-static de- 
pinning exists, with velocity 



f-fc 



(3.17) 



k(7o + Vo) - Vo ' 
until the critical value of the interlayer coupling is reached, 

(3.18) 



Vo = Vc 



1 - K 



7o 



Here the reentrant (or hysteretic) branch appears. This is qual- 
itatively similar to the mean-field picture. One can relate for- 
mally (1 — k) — » M(c) ~ h 2 /c 2 which for small disorder has 
the same form as d3.161 l. if the elastic coefficient (set toone is 
this Section) is restored and one identifies A"(0) — ► h 2 . An 
interesting question is the nature of the elastic to hysteretic 
transition at ?; c . Expanding ( 13.13t in powers of v yields the 
equation of motion near the critical point: 



(Vc -V)v = f- fc + 2^(7n + 7i2)A'"(0+) 



1 



+ 0(v 3 ) 



(3.19) 

As one can see on Figure 13.11 the transition is continuous if 
A"'(0 + ) < and v ~ VT~ r Tc at the transition 3 . Such a sce- 
nario may hold for the non-periodic, random manifold class 4 . 
A series of higher multicritical points should exist, associated 
to correlators with leading behaviour A'" +1 ^ (0 + )v n . For the 
periodic scalloped potential n = 00 and the transition exhibits 
a jump, or a quasi-jump (inverse logarithm) as in mean-field, 
illustrated in Fig |3.2l 



3 From the factor J q 6 one could identify d = 6 as a critical dimension 
for the tricritical point, and find that the terms Dud^u and Bu(dtu) 2 
both become relevant there. However one should remember that A'" (0) is 
irrelevant in d = 6. Whether this modifies the exponents and leads to new 
universality class is left for future study 

4 For d > 4, if th e bar e disorder is strong enough, A (it) develops a cusp, 
see Appendix in ll03ll . 
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(ii) d < 4: the friction coefficient of the center of mass 
decreases with scale as: 



1+ (L)^( lQ+m )(L c /L) 2 - 
It reaches values near zero at a scale 



Lp\ — L c 



7o + Va 



(3.20) 



(3.21) 



which diverges as 770 — * 0, and which we term the "plastic 
length". Thus the depinning transition of a system of size 
L < L p i remains similar to the standard (finite-size) elas- 
tic depinning of a single layer, while systems with L > L p \ 
cannot be described by single-layer elastic depinning. It is 
then likely that the system breaks into domains which can de- 
pin and move independently. The full collective dynamics at 
scales L > L p \ however remains to be understood. This insta- 
bility of the elastic depinning at finite scale is an effect beyond 
mean field. 

Another important length scale is associated to a non-zero 
velocity. For single-layer elastic depinning it is 



AM 7 
7o v 



(3.22) 



It is such that vt — Tf, where r is the time scale diverg- 
ing at depinning and rj the correlation length of the disorder, 
equal to the period a (here set to unity) for the simplest CDW 
class. Beyond that scale the effect of quenched disorder is 
washed out into an effective thermal noise and the motion is 
uncorrected. Equating the two scales L„(7o) = L p \ defines a 
characteristic velocity scale: 



A 2 



'7 



7o + Vo 



(3.23) 



below which plastic effects cannot be neglected. The be- 
haviour of the system at and beyond that scale still needs to 
be elucidated. 

Since we found that the FRG fixed point of elastic depin- 
ning is dynamically unstable (to one loop) to the viscous cou- 
pling we now investigate the phase diagram of the moving 
phase. 

B. Functional RG 

Let us now derive and analyze the FRG equations to 1-loop 
order at non-zero velocity. For pedagogical purposes, we use 
a Wilson scheme i.e. we compute T[u, u] to one loop using a 
cutoff A; = Ae~ l and write RG equations as the cutoff is var- 
ied (i.e. integrating over a shell using J Aie - M<g<Al /(?) = 

S^Af /(A;)d0- A method which can handle higher loops, 
based on a non-zero mass scheme, is presented in the next 
section. Here we restrict to the periodic problem and choose 
units such that the period is one. 

The standard single-layer result for the correction to disor- 
der upon integration over the shell can be expressed as 



d t A(u) 



eA(u) - \ 



A(u) - A(0) 



(3.24) 




FIG. 3.1: schematic v — f curve corresponding to Eq. | |3. 19t 
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FIG. 3.2: schematic v — f curve corresponding to a sharp transition 



where one has defined A(u) — Sd^ e A(u). This result holds 
in the limit of zero velocity v = + . As is well known, it 
results in a non-analytic correlator beyond the Larkin length 
L c . We note that a non-analytic A(u) decreases 7, while an 
analytic correlator would increase it. We denote 



<7 = A"(0+) , a = A"(0+) 
The family of quadratic correlators: 

'1 



A(u) 



6 



it(l — u) 



(3.25) 



(3.26) 



for < u < 1, periodically continued to all u, is preserved 
by the FRG flow, with dio = ea — 3<r 2 . It is realized by 
a scalloped potential, or more generally by uncorrelated pe- 
riodic shocks, and contains the universal fixed point of the 
random periodic (RP) class: For e = 4 — d > it flows to 
a* = 2 — z c \ = e/3. For d = 4 the fixed point is at zero but 
the slow asymptotic decay a ~ 1/(31) is universal. 

Inserting formula d3.6t for the response function and the 
Fourier series A(u) = J2 P e l2vpu A p (over integer p) one 
finds the correction: 



£711 





1 






Jq 


_2v 2 > 




-I 


du 







1 



A"(u) 

A"(o+y 



2u 2 7 



(3.27) 

(2tt P ) 2 A p 
(x — i2np) 2 a 
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For the scalloped potential family ( 13.26b A p = (1 

5 p o)cr/ (2ttp) 2 , and $(.t) reads 



*(a) = --a 



1 



[2 sinh(x/2)]' 



(3.28) 



In the sequel, we use the scalloped family ( 13.26b and the form 
d3.281 l. From d3.27l i and ( 13.28b one obtains the RG equations 



9i7n = -o"7ii + 7rr 



1 



A, 4 



7 + sinh 2 (^_) 



7 sinh 2 (^) 
(3.29) 

and of course <9;7i2 = 0. For v = + it reproduces the elastic 
depinning RG equation 3; 7n := (z c \ — 2)jn = — tryn which 
at the fixed point yields the dynamical exponent z c \ = 2 — e/3. 

One can see from formula ( 13.29b that a non-zero velocity 
v > tends to cut the flow of 711. This is a usual effect in 
the case of elastic depinning (770 = 0) associated, in that case, 
to the single length scale L v (j ) defined in (13. 22b . Here there 
are a priori two length scales, associated to the two modes u + 
and u~ . The effect of disorder is washed out only for scales 
larger than both lengths, i.e. if 

7±ve 2 7A 2 > 1 , L> Lf . (3.30) 

Then the equation reduces to di 711 = ~^> 7+ ^ v i &~ 41 - The 
difficulty is that Lf are not simply equal to L v (7° ) since the 
7±(0 do not behave as the single-layer coupling (it does only 
for scales L < L p \). In fact, 7+ (7) may vanish at some scale, 
hence the condition j+^ve 21 /A 2 ^> 1 may never be fulfilled, 
at any scale. A more careful analysis, performed below, is thus 
required. 

One notes that (13. 29b is the derivative 9;7n = 9(7+ = 
—d v dif with: 



dif = 7ii.5"!; 



aA 2 



A 2 

coth(- 1 



A 2 

coth( 7 ' 



(3.31) 

from which the velocity-force characteristics is obtained as: 



dld t f 



In the limit v = + one recovers dif = — ^<5A 2 which in- 
tegrates to — ff' sl = — jctA 2 , the critical force of a single 
elastic layer. One notes the general relation, 



d v f(v) = 7+ 



(3.33) 



valid for I = 00, or for any intermediate scale, if one defines 
a finite-scale curve for f{v) by setting the upper integration 
bound to I in ( 13. 32b . 

We now study the flow of 711 which depends on 7+ = 
7n — Vo an d 7- = 7+ + 2^o- We recall that the starting value 
i s 7ii = 7o + Vo- If me velocity is large enough, although 7+ 
decreases upon renormalization, the corrections may be weak 
enough so that it remains positive, even for d < 4. In the lat- 
ter case, there should always be a critical velocity v c such that 
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FIG. 3.3: v(F) curve obtained by integration of the one-loop equa- 
tions in the text, for A = 1, a = 1/3 (e = 1) 770 = 10, 70 = 1 



7 + (7 = 00) = 0. For v > v c the v(f) curve is well-defined 
and continuous. For v < v c there is no moving solution such 
that v(f) has a positive slope. Hence in a fixed-applied-force 
ensemble there is a jump to the pinned phase as the force is 
decreased. At v = from ( 13.331 ) the slope of the v(f) 
curve is infinite. This corresponds to the minimal force /< at 
which the jump must occurs. By contrast, when the force is 
increased in the pinned phase the critical force is /> = f^' al . 
It corresponds to the maximal force at which the jump upward 
in velocity to the moving state must occur 5 . At v = w+ the 
length L+ is infinite. This suggests that motion should be cor- 
related on all scales, and that this point is very much like a 
critical point where scale invariance holds. An example of a 
v(f) curve predicted by the one-loop FRG is given in Fig. 13.31 
To estimate the jump velocity v c it is simpler to first study 
a model where the bare value of 7+, 70, is already small com- 
pared to 7_, i.e. 70 <C i]o- Then equation ( |3.29b can be ap- 
proximated by: 



<9z7+ = -5ryo 



16r)ov 2 sinh ( 



2/ A 2 e- 



(3.34) 



4)70 v 



(3.32) This is integrated into: 



1 / A 2 
7 + (Z = 00) = 70 - -arj a H[ 

2 \^i]ov 



Jo V 



1 - 



sinh (y) 
1 + x coth(a;) + In 



(3.35) 



(3.36) 



^sinh(x) 

with H(x) w ln(2x) — 1 at large x and H (x) ~ ^- at small 



5 In some cases it was observed that the jump can occur before these extremal 
values, either due to finite-size effects or due to a dynamical instability, 
which is beyond the present description. 
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x. The critical velocity v c is hence determined by 

\4r) v c J 0-770 
which gives the asymptotic behaviour: 



(3.37) 



A 2 

VQVc 

A 2 



i exp 



370 



, ^ « 1 (3.38) 



270 

°Vo 
1/2 



7o 



7o 



> 1 



(3.39) 



Here we have assumed a to be scale independent, hence a 
reasonable value for it 6 is a = a* = 2—z e i = e/3+0(e 2 ). At 
the fixed point in d = 3, 2, 1 the second regime is the relevant 
one and gives the value of the critical velocity for large rjo /70. 

To estimate the critical velocity when i]q / 70 is small, one 
must first integrate the flow up to scale l\ at which 7+(?i) = 
krjQ = fc/(fe+2)7_ and k a number smaller than unity. Within 
this scale we can approximate 711 (0 = (rjo + 7o)e _cr ', which 
yields (?7o+7o) e ~ <T ' 1 = (k+l)r)o. The length scale L c e l± is of 
the order of the plastic length L p \ introduced above. Beyond 
that scale one can apply the previous analysis 



1 / A 2 



which yields the estimate 



rjoVc 
A 2 



(fc + 1)7/0 
7o + Vo 



2/ct 



m- 1 {2k/d) ■ 



(3.40) 



(3.41) 



Hence we find that the critical velocity vanishes as ijov c ~ 
(%/7o) 2 ^ 2_z ' 1 m the limit of small viscous coupling, consis- 
tent with the estimate ( 13.231 ) for the scale at which plastic ef- 
fects become important. The present 1-loop analysis indicates 
however that the jump is always non-zero 1 . 

It is instructive to compare with the predictions from mean- 
field theory (MFT) recalled in Section UTAl In the regime of 
large viscous coupling rjo ^> 70, one sees that formula d3.39t 
is very similar to the mean-field prediction 



VoV c 



(c + h ) 



Vo 
37o 



1/2 



(3.42) 



if one identifies c + h — > Hence the 1-loop FRG 

result, taken in the limit of large 770, is very similar to mean- 
field theory (MFT) even for d < 4, with the difference that 



1 Given the assumption 70 -C »?o the first regime exists only for small a 
which is either e — » 0, or if bare disorder is very small until the scale 
which controls the jump. 

One notes that the flow of the disorder correlator, which is too complicated 
to analyze here, is also cut by velocity at the scale max(L„ , L^)- Hence 
above that scale the parameter a cannot be assumed to take its fixed-point 
value and instead will decrease to zero. Since the effects computed here 
occur below these scales, one expects at most a change in the prefactors 
from taking these effects into account. 



the disorder parameter flows to a universal fixed value a* . In 
the other limit of small ratio 770/70, the result is very different 
from MFT because of the strong renormalisation of the in- 
layer friction coefficient, and the threshold i] c which exists in 
mean field is zero for d < 4. 

It is also instructive to study the FRG flow for d = 4 and 
d > 4. For d = 4 + e and a scalloped potential one has 
a = <xoe _d , hence one finds at zero velocity 711 (/) = (70 + 
770) cxp (-^(1 - e~ d )) and 7+(0 = 7n(0 - Vo- There is 
thus a threshold for the jump in the v(f) curve; it occurs only 
for 770 > ij c with 



Vc 



To 



1 



(3.43) 



?7 C becomes very small as d — > 4 + . For d = 4 one has a = 
1/(3?), hence 711 = (70 + Jyo)/? 1 ^ 3 and there is no threshold, 
r] c = 0. The plastic length scale however diverges extremely 
fast as L p \ = L c cxp((7o/?7o) 3 ) for small t)q. 

The analysis of this section used that j 12 is not corrected. 
We now turn to a two-loop analysis to check whether this 
holds to higher orders. 

IV. ANALYSIS INCLUDING 2-LOOP 
CORRECTIONS 

In this Section we compute the corrections which arise at two 
loop around the quasi-static elastic depinning transition of the 
single layer at / = /® • . The calculation is performed in 
the limit v — > + . The FRG flow is discussed separately for 
the non-periodic and for the periodic cases. Possible conse- 
quences at non-zero v are discussed in each section. 

The natural setting for higher-loop calculations is to use a 
mass as an infrared cutoff. It amounts to adding the force 
vector m 2 (w{t) — u % x t ) to the r.h.s of the equation of motion 
dl.5l l. It describes two layers both pulled by a spring attached 
to a point at position w(t) which perfor ms qu asi-static for- 
ward motion. In that setting, it was shown Il04ll that the force 
correlator A(u) computed in the FRG is an observable related 
to the mean-square center-of-mass fluctuation around w(t) in 
each layer. One introduces the rescaled correlator 



A(V) = C d m" e+2C A(7im- c ) 



(4.1) 



where C d = el 2 = e j k (k 2 + l) 2 for e = 4 - d > 0. One 
finds that A(w) converges to a fixed point, and to 1-loop order 
it reproduces the Wilson approach. 

A. 2-loop FRG equations 

The 2-loop FRG flow-equation for the disorder is taken to be 
the same as the one derived in ifTcill at the quasi-static depin- 
ning transition: 

d e A(u) = (e - 2C)A(u) + C«A'(u) 
1 

~2 
1 

+ 2 



(a(u) - A(0i 
(A(«)-A(0)) A'(uf 
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-iA'(0+) 2 A'» 



(4.2) 



where dg := —md m . As explained there, the derivation of 
this FRG equation at two loop (especially the last term) relies 
on the Middleton theorem 112311 which states that if all local 
velocities are positive at some time, they remains so at all 
times. In the two-layer viscous problem this property does 
not hold stritly, as backward motion of one layer is sometimes 
observed. The present calculation hence assumes that these 
effects can be neglected at large scale near the quasi-static de- 
pinning, and to this order. 

The corrections to the friction coefficients are computed in 
Appendix[B] They read: 



_ 7l2 A'(0+)A'"(0+) 
o e 712 = — - log 



7n + 7i2 



7n — 7i2 



dg 7n = 7n 



(4.3) 

-A"(0) + A"(0) 2 + A"'(0)A'(0) (§ - In 2) 



37i 2 A'(0+)A"'(0+) 



log 



7n +7i2 



+7nA'(0 + )A'"(0 + )log 



7n — 712 

_^ 

7i 2 i 



(4.4) 



The calculation was performed in the physical domain 
7+,7- > 0. For mainly illustrative purpose, an analytical 
continuation was performed to the domain with negative fric- 
tion coefficients, which yield the absolute values above. We 
find however that whenever the coefficient of the log terms are 
non-zero, the solution of the flow, obtained below, remains in 
the physical region. 

It turns out that the two combinations A"(0 + ) and 
A'(0 + )A"'(0 + ) which appear in these equations are univer- 
sal numbers which can be related to the roughness exponent, 
£ (using derivatives of ( 14.2b at u = 0), independently of the 
precise form of the fixed point: 



A"(0) 



A'(0)A( 3 >(0) = 
Here we have defined 



1-Ci C[_3Ci-3C 2 + 2^ 2 
3 6 

(l-Cl)Cl 2 



9 



e 2 + 0(e 3 ) 
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C = Ci£ + C 2 e 2 + 0(e 3 ) 



(4.5) 



(4.6) 



B. Non-periodic problem 

As was shown in fTill for a wide range of microscopic dis- 
orders, there is a unique elastic-depinning fixed-point, calcu- 
lated there, and identical at one loop to the random-field (RF) 
disorder class. At this fixed point one finds £i = 1/3 and 
C 2 = 1/(27V2 7 ) with 7 = 0.5482228... This yields: 

dan = (-0.0432087e 2 - ^ 7n 

+ i.o g (i- 1) ^ 




FIG. 4.1: Flow of 711 and 712, for e = 1, as a function of 711 (x- 
axis) and 712 (y-axis). The separatrix is the diagonal line (orange) 
7+ = 711 + 712 = 0. All physical initial conditions , corresponding 
to 7+ > 0, remain physical. The shaded pink region corresponds to 
unphysical initial conditions 7+ < 0. 



<9<?7i 2 



108 



e log 



7n + 7i2 
7n - 712 
7n + 7i2 
7ii — 712 



712 e 



712 



(4.7) 
(4.8) 



We integrated the flow-equations numerically. The result is 
given in Fig. 14.11 for e = 1 and in Fig. 14.21 for e = 4, to 
illustrate how the flow changes with e. Looking carefully, one 
sees that starting in the physical region 7+ > 0, the unphysical 
region 7+ < is avoided. One also sees that 7+ approaches 
zero quickly, at least for small e. We now confirm these two 
findings analytically. To do so, we change variables to 7+ = 
7n +712 and 7_ = 711 — 712, see Eq. ( II.6I 1. We are interested 
in 7+ w 0. There the flow-equations become 



da- = - 



[ ^ij e- 0.008768e- jjj ] 



0.008768e e-y- 



The second equation has the solution 

-J-(l) = c -( 1 /9+0-008768e)rf 7 _(-Q^ _ 



(4.9) 



(4.10) 



The solution for 7+ is easiest expressed as a function of 7_, 
instead of I: 



d 7+ (7-) M 7 ^) 

d 7 _ 12 + 0.94697e 



The ratio 



r(7_) 



7+ (7-) 



(4.11) 



(4.12) 
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FIG. 4.2: Flow of 711 and 712, for e = 4, as a function of 711 
(x-axis) and 712 (y-axis). The separatrix (orange) is the line 7+ = 
711 + 712 = 0. The pink region corresponds to unphysical initial 
conditions 7+ < 0. 



satisfies a closed flow equation as a function of 7_ : 



dr(7_) 
din 7_ 



r(7_) - 1 



elogr(7_) 
12 + 0.94697e 



(4.13) 



For all relevant values of e (0 < e < 4), there are two solu- 
tions: r = 1 (unstable) and a non-trivial (r* <C 1) solution 
of 



1 = 



elogr* 



which yields: 



exp 



12 + 0.94697e 



12 + 0.94697e^ 



(4.14) 



(4.15) 



The eigenvalue of the flow close to r* is at leading order 

£ 2 e 12/ e 



V 



108 



(4.16) 



Thus for e small, this fixed point is very attractive. This is the 
fixed point obtained numerically above. It has the property 
that 7 + remains strictly positive. 

From d4.10| ) we extract the dynamical exponent associated 
with 7_ : 



Aplastic 



2-1 
9 



0.008768e^ 



(4.17) 



Since 7+ ~ r*7~, it has the same dynamical scaling, and the 
above Aplastic is indeed the critical exponent for the dynamics 
of both modes. 



Hence within the two-loop analysis, and the stated assump- 
tions, one finds a fixed point for the case of non-periodic disor- 
der. The dynamical exponent at this new fixed point deviates 
even at leading order in e from the standard elastic depinning 
value: 



^clastii 



2e 
~9 



0.0432087e J 



(4.18) 



Compared to one loop, the two-loop corrections appear sin- 
gular, as seen from the ln(7_/7 + ) factors in the corrections 
to friction. As a result their magnitude is drastically enhanced 
above the plastic length L p \ from an expected 0(e 2 ) to an ac- 
tual 0(e). The term A'(0 + )A"'(0 + ) 111(7-/7+) in the correc- 
tion to 712 in Eg. (14.3b is in effect replaced, upon integration 
of the flow, by A"(0). This results in a value for 2 — z twice 
smaller, to leading order, than the usual elastic fixed point. 

To summarize, the 1-loop analysis showed that 7+ becomes 
very small near the plastic length, and provided a scenario for 
scales larger than L p \ which could sustain only a moving state 
at v > v c . Although we did not perform the analysis for the 
non-periodic case in detail we do not expect a difference at one 
loop. The present analysis - in the non-periodic case - shows 
that additional physics occurs at two loop. It suggests that a 
v = + state may still be possible. From the above analysis, 
one could surmise that it results in a very abrupt, almost ver- 
tical v(f) curve (since 7+ is found to converge very rapidly 
to a very small value) which is not strictly a jump, although 
it may look like one in a numerical calculation or an experi- 
ment. This "quasi-jump" would occur near the critical force 
of the elastic system, at variance with one loop. To confirm or 
infirm this scenario one would need to include the effects of a 
non-zero v, and a possible violation of the Middleton theorem 
within the two-loop theory, a challenge left for future work. 

It also remains to be investigated to which extent the present 
analysis can be trusted in the region where 7+ becomes very 
small, i.e. the region where ln(7_/7+) becomes of order 1/e. 
We expect that in that region terms such as u 2 in the equation 
of motion may become important. Such effects are presum- 
ably correctly resummed in the two-loop corrections and may 
explain why 7+ remains positive. However since the count- 
ing of order in e becomes unconventional if one follows the 
flow further in that region, there is no guarantee that higher 
loops may not lead to even more singular terms. In the best- 
case scenario only the 0(e 2 ) term in ( 14. 1 8b would be changed 
by higher-loop corrections. Although the present results hint 
at a new depinning universality class with a dynamical ex- 
ponent z = z p \, a deeper understanding of the behaviour of 
the system in the plastic region seems necessary before a firm 
conclusion can be drawn. 



C. Periodic problem 

The case of periodic disorder is also challenging. The quasi- 
static depinning fixed point has the form A(u) ~ u(l — u), as 
in d3T26b with 



e 

18 



(4.19) 
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and is expected to maintain that form to any order in e. If the 
system is exactly at its fixed point, then, since A"'(0 + ) = at 
this fixed point, the flow-equations for the 7's read 

(4.20) 

- 7ll (2-z el ) (4.21) 



(4.22) 



Hence there are no drastic effects of the two-loop corrections, 
apart from changing the value of z: 7ll always decreases as 
in 1-loop, 7+ vanishes at some scale, and the 1-loop analysis 
remains at least qualitatively correct. Hence this confirms the 
1-loop approach. 

It is less obvious to understand the situation where the sys- 
tem is not exactly at its fixed point, but converges to it, i.e. 
A"'(0 + ) ~ e~ el . Inserting this behaviour in the above two- 
loop equations still results in drastic effects, i.e. 7+ never 
crossing zero, again due to the logarithmic divergence of the 
corrections in that region, as for a non-periodic problem. The 
discontinuous behaviour between a zero and a small non- 
zero A"'(0 + ) remains to be understood. One scenario which 
would save the agreement with the 1-loop approach is that 
other irrelevant operators than A"'(0 + ), neglected in the two- 
loop treatment of the periodic class, are equally important and 
modify the result back to (|4.2Q| i. More work is clearly needed 
to settle these issues. 

V. TOY MODELS WITH 2 PARTICLES 

To gain insight on some of the issues arising in the dynamics 
of coupled elastic layers it is instructive to study the model 
in d = i.e. a toy model with two particles. This appro ach 
has proved useful for the elastic-depinning problem Il04ll . in 
particular in clarifying the information contained in the FRG 
functions. As we show below, a variety of behaviors arises 
already for two viscously coupled particle. Here we focus 
on the simplest situation of two particles in a periodic one- 
dimensional landscape driven by a force, and leave for future 
work the interesting non-periodic case, as well as driving by 
a spring (which is more suitable for comparison with FRG). 
The model is thus the d = version of ( 12. It , with a pinning 
force h l Y(u l — (3 l ); we choose h 1 = h 2 for simplicity. The 
random phase can be eliminated by a shift of the u l , hence 
it is sufficient to study the case of two particles in the same 
landscape (up to a change in initial conditions). We first study 
smooth disorder, and then a scalloped landscape 8 . 

A. Smooth potentials 

We now study the following model: 

7«i = r]{u2 - ill) + f + 4>(ui) 



7«2 = T)(U1 ~ U2) + f + <i>{u2) 



(5.2) 



' If an additional self-consistency condition is imposed, these models can 
also be used to implement a third mean-field approach, discussed at the 
end of Section lHAl 



In this Section we adopt slightly different notations for center- 
of-mass and difference coordinates: 



y 



Ul + u 2 



(5.3) 
(5.4) 



X = Ul — U-2 ■ 

In these coordinates, the equation of motion becomes: 

y = F-±[<j>(y + x/2) + <f>(y-x/2)] (5.5) 
x = a [<t>(y - x/2) - cf>{y + x/2)} , 

where we have defined: 

7 

a = ;r~ 

7 + 2ry 



F = 



(5.6) 

(5.7) 

(5.8) 
(5.9) 



For definiteness we consider the family of periodic-force land- 
scapes: 



(«) = 



Pi sin(27ru) + P2 sin(47ru) 



- m% +p 2 2 + (i + 125?) V( 1 + 32 ^) 

(5.10) 

They are normalized such that if one takes \pi\ = 1 (the 
standard choice made in the following) the single-particle 
critical depinning force is F^ p = ±1 (i.e. max(4>(u)) = 
m&x(—(j)(u)) = 1) for any p^. It turns out that the single- 
harmonic case is non-generic and one needs to include at least 
one other harmonics, i.e. p 2 7^ 0. 

We have integrated these equations numerically and plotted 
9 the resulting flow in Figs. 15.1115.21 and 15. 31 for various val- 
ues of / and p2- The center-of-mass coordinate y is plotted 
along the vertical axis, the relative displacement x along the 
horizontal axis. 

It is instructive to start with the case 77 = (a = 1) of two 
uncoupled particles, given in Fig. 15.11 The vertical trajecto- 
ries along the y axis at x — or x = 1 correspond to the two 
particles either in the same position or shifted by one period. 
As the force is decreased below threshold (right to left) a pair 
of fully-attractive and fully-repulsive fixed points appears on 
these axis. The total phase space for F < F c sp = 1 is frag- 
mented in pinned regions which flow to one of these "pinned- 
phase" fixed points, corresponding either to the two particles 
pinned in the same well or pinned in two wells shifted by one 
period, depending on the initial condition. Note also the other 
zero-force fixed point which has one attractive and one re- 
pulsive direction and corresponds to one particle in a stable 



9 We are grateful to Alan Middleton for clarifying remarks during the analy- 
sis of these flows. 
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FIG. 5.1: Uncoupled particles below depinning (F = 0.85, left), at depinning (F = 1, middle) and above depinning (_F = 1.2, right); a = 1, 
pi = 1, p2 = 0. We always plot x to the right and y to the top. Separatrices for the different attractive regions below threshold are drawn in 
green. 




FIG. 5.2: Viscously coupled particles (a = 0.2) below depinning of the uncoupled particles (F = 0.6). The anharmonic coefficient p2 differs 
from left to right: p2 = (left); p2 = 0.5 (middle) and p2 = —0.5 (right), pi = 1. We plot 20 sample trajectories starting from y = 0, and 
equally spaced in x. One sees that for p2 > more trajectories converge towards the unique stable solution (yellow). In the case of p2 = 0, 
there is a family of periodic solutions, of which we have plotted three. 



equilibrium position at the bottom of one well and the other in 
an unstable equilibrium position at a hill top. This fixed point 
controls the separatrices of the flow. This structure, obvious 
in the absence of a coupling, will persist, with some modifica- 
tion, for non-zero rj. 

Interesting physics happens when the viscous coupling rj is 
increased. The case a = 0.2 is shown in Fig. 15.21 Exactly 
along the axis x = and x = 1 the equation of motion has 
not changed and the same attractive pinning fixed points are 
present for F < F* p = 1. However, unbounded motion is 
now possible for smaller forces F c < F < F c sp = 1, and 
takes place away from the axis. The force chosen in Fig. 15.21 
is F = 0.6. On the left figure the case p2 = is repre- 
sented. One can easily see that it is fully integrable and that 
each trajectory in the central region is exactly periodic and 
crosses the y = 0, 1 axis at the same x. The region where this 
flow occurs is delimited by the separatrices which meet at the 



above-mentioned zero-force saddle points. Hence one sees 
clearly that the phase space splits into a pinned region and a 
flowing periodic region. In the case of a pure sine (p2 = 0), 
this region is made of an infinity of neutral periodic trajecto- 
ries (with zero Lyapunov exponent). In the more generic case 
P2 7^ 0, the flowing region contains a single periodic trajec- 
tory. This trajectory is either attractive (case p2 = 0.5, figure 
15.21 in the middle) or repulsive (right part of figure 15.21 with 
P2 = —0.5). It is easy to prove from the symmetry properties 
of the flow that the Lyapunov exponent is reversed when the 
sign of the force landscape is reversed <f>{u) ~ > —<fi(u) 10 . In 



10 Denoting ■u I (t, /, pi,p2) the solution of the equation of motion 
- for some given but unspecified initial condition, one sees that 



«**(*, ~f,Pl,P2) 



-U t (.t,-f,p 1 ,p 2 ), U z (t,f,p!,p2) 



1/2 
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FIG. 5.3: Viscously coupled particles (a = 0.2) above (F = 0.7, right), at (F = F c = 0.522264, middle) and below depinning (F = 0.35, 
left). The anharmonic coefficient p2 — 0.5. We plot 20 sample trajectories starting from y = 0, and equally spaced in x. For F = 0.7 (right), 
we plot the unique stable solution (yellow). Even for this large F, one sees the convergence to this stable solution. 



the repulsive case, any particle in the region apparently flow- 
ing on the figure eventually gets pinned at some larger y, after 
visiting a few cells; the basin of attraction of the flowing phase 
has measure zero. This is an example where a non-trivial pe- 
riodic stationary state exists, but is dynamically unstable. On 
the contrary, in the case p2 = 0.5 (middle of figure [V~AT > a 
distinct flowing phase exists, and its properties are dominated 
by a unique attractive periodic trajectory, and e.g. the average 
velocity is given by the inverse period of this trajectory. 

Finally Fig. |5.3| illustrates how the periodic orbit in the mid- 
dle, hence the moving phase, disappears when the force is re- 
duced below F c = 0.522265, leaving only a pinned phase for 
F <F C . 

We can now analyze the resulting v(f) curve. The v(f) 
curve for the pure-sine model is indicated schematically on 
the left of Fig. 15.61 and is non-generic, as discussed above. 
In the case p2 ^ there are two branches corresponding to 
the two steady states, one (labelled 1) corresponding to the 
trajectory along the x = 0, 1 axis, i.e. the single particle v(f) 
curve, and the second (labelled 2) corresponds to the periodic 
orbit near the middle of the figures , which generally has a 
higher v(f) curve. If the second is repulsive (p2 = —0.5), 
then the trajectory along the x = 0, 1 axis is attractive: the 
global v(f) curve then coincides with the single-particle one 
and there is no hysteresis (middle of Fig. 15.6b . If the second is 
attractive (p2 = 0.5), the Lyapunov exponent of the periodic 
trajectories are inverted 11 and the global v(f) curve follows 



~Pi,P2), and u % (~t,f, -pi, -p 2 ) = -u l (-t, f,p!,p 2 ) = 
1/2 — t, f, —p\,p2). This last property implies that the two leftmost 
figures in Fig. IV A| can be deduced by symmetry and that the Lyapunov 
exponent on the periodic trajectories (which are globally preserved by the 
symmetry) are reversed in sign. The pinned fixed points however remain 
attractive and are simply exchanged by this symmetry (they are not indi- 
vidually preserved). 

Note that while the line x = 0, 1 are always attractive in the vicinity of 
the pinned fixed points for F < 1, it becomes - in that case with pi = 



the second branch. In that case there is a hysteresis as the force 
is varied adiabatically; this is shown in the right figure. Upon 
decreasing the force from a large value the system follows 
the attractive trajectory in the middle until it disappears at F c 
and the velocity vanishes. But if the force is increased from 
a value smaller than F c , it can be seen from the left plot on 
Fig. |5.3| that it first converges to a pinned fixed point along the 
axis x = 0, 1. Since these fixed points remain attractive up to 
F = F c sp = 1, the velocity remains zero until that force and 
then jumps to the stable moving state. 

The question of whether a jump exists in the descending 
curve can be settled by analyzing how the periodic trajectory 
disappears at F = F c . It can be seen from the middle plot on 
Fig. 15.31 that this occurs abruptly, but that the period diverges 
at F = F c + as the system spends more and more time near the 
zero-force saddle points. These hence play an important role 
in the transition at F = F c . A simple argument indicates that 
the time spent near these points increases logarithmically, as is 
verified by the numerical integration of the flow in Figs. 15.41 
15.51 Hence, although this system exhibits hysteresis in the 
case p 2 = 0.5 it does not exhibit a velocity jump along the 
descending branch. Note that the critical behaviour at F c is 
different from the single-particle case v ~ (F - i^P) 1 / 2 , due 
to the zero-force saddle-point mechanism. 

B. Scalloped potential 

Here we consider the two-particle toy model for a piecewise 
parabolic (scalloped) potential, corresponding to a piecewise 
linear pinning force with jump discontinuities at the bound- 
aries of each period. The equations of motion for the position 



0.5 - repulsive for F > 1 when the flow starts along this line. This is 
again a consequence of the symmetry properties mentioned above which 
inverts the Lyapunov along a periodic trajectory globally preserved by the 
symmetry. 
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FIG. 5.4: Velocity uasa function of F. The parameters used are 
Pi = 1,P2 = 0.5, a = 0.2. 
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dinates, the equations of motion are 

1 n + m 

iv = -y + f + 2 + — 2 — 



jx = —ax + a(n — m) 



(5.14) 
(5.15) 



for n±m < 



y < n±m + 1 and n - (m + 1) < x < 1 + n — m, 
with a given in Eq. (|5.7| i. Our goal is to identify the stable 
periodic orbits for this model and calculate the corresponding 
period or its inverse, the mean velocity. The regions in coordi- 
nate space corresponding to the various periods of the pinning 
potential are shown in Fig. 15.71 Consider a particle that starts 
at point A with [x(0), y(0)] = (xq, — xq/2) in the region of 
the pinning potential corresponding to (n,m) = (0,0). The 
boundaries of this region are defined by — x/2 < y < l+x/2, 
for -1 < x < 0, and x/2 < y < 1 - x/2, for < x < 1. 
This particle will travel across the (0, 0) region to point B in 
a time t\ (xq) and then across the (0, 1) region to a point C in 
a time t2(xo), according to 



A 



Xq\ t!(x ) 



D 



x(h),l 



\x(h)\ 



C= (x{ti + 1 2 ), 1 + 



|a>(ti + ta)| 



(5.16) 



The case of a scalloped pinning potential can be studied ana- 
lytically since the equations of motion are linear within each 
pinning period, with jump discontinuities in the velocity at 
boundaries of the pinning regions shown in Fig. 15.71 



FIG. 5.5: Velocity v as a function of -l/ln(F - F c ), with F c = 

0. 522265. The fit-function is v = 0.0227714 - 0.718327/ ln(F - 
F c ). The linear fit is excellent. A (much worse) fit to a power-law 
would give an exponent of about 0.1. The parameters used are pi — 

1, p 2 = 0.5, a = 0.2. 



a. Periodic orbits We wish to determine the values of xq 
that correspond to periodic orbits as defined by the fixed point 



x(ti + t 2 ) = x'(x ) = x 



(5.17) 



The period of such orbits is t\ + t 2 and v = l/(ti + 1 2 ). It is 
convenient to introduce a new notation: 



of the particles are 



r](u 2 - ii! ) + /+-+ n - ui (5.11) 



1U\ 

JU 2 = - u 2 



-+m-u 2 (5.12) 



for n < ui < n + 1 and m < u 2 < m + 1, with n and m 
integers. 

When rj = 0, the particles are decoupled, and the dynamics 
can be determined exactly. Each particle is pinned for / < 
1 /2. For / > 1/2 there is a unique periodic orbit of period 



v 1 V/-1/2 



(5.13) 



that diverges linearly as / — >■ (l/2) + . No periodic orbits exist 
for / < 1/2 and the system does not exhibit hysteresis. 

To consider the case 77 7^ 0, we introduce center-of-mass 
and difference coordinates as in Eq. ( 15.3b . In these new coor- 



z 1 (x ) 
z 2 (x ) 



-tl/7 
-t 2 /7 



with 



7ln(l/ziz 2 ) 



(5.18) 
(5.19) 



(5.20) 



The dynamics from (xq, |a;o|/2) to (x 1 , 1 + |a;'|/2) can be ex- 
amined analytically since the equations of motion are piece- 
wise linear. It is determined by 



fz2-\z a 2 + ] ^z a 1 {z 2 + z a 2 ) = f- 1 - 



with 



xqzIzI 



F0| 



(5.21) 
(5.22) 

(5.23) 
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FIG. 5.6: Schematic v(f) curves corresponding to the three cases discussed in the text: (i) left: single sine force landscape (ii) middle: the 
non-trivial periodic orbit 2 is repulsive and the v(f) curve is the same as for a single particle 1 (iii): left: the periodic orbit 2 is attractive and 
the v(f) curve exhibits hysteresis as discussed in the text 




FIG. 5.7: The figure shows the boundaries of the regions where the 
relative and center-of-mass velocities of two particles in a periodic 
scalloped potential have jumps. The horizontal and vertical coordi- 
nates are the relative and center-of-mass position of the two particles, 
respectively, as defined in Eqs. J5,3b . The equations for the straight 
lines bounding the region corresponding to n = 0, m = are indi- 
cated in bold. 



We now look for a periodic solution or fixed point as defined 
by Eq. (5). Then Eq. ( 15.23b gives (provided uq ^ 0) 



Fol = 



l-*2 
1 - z\z% 



(5.24) 



Substituting this in Eqs. (15.21b and (15.22b we obtain 



z 1 (f + l/2)- { \;Z ){z l + J ) =f-l/2, (5.25) 



z 2 (/ + l/2)- 



2(1 - z-z-) 

(l-sf)(z 2 + *g) 
2(1 - z^) 



f-1/2. (5.26) 



These two equations are symmetric in z\ and z<z, indicating 
that the solution must satisfy z\ = z-2 = z. There is a fixed 
point x* of Xq where the system undergoes a periodic orbit of 
period 1/v = —2 ln(z), with 



1 

1 + z a 



(5.27) 



1.0 r 




FIG. 5.8: The velocity-force characteristic obtained by inverting 
Eq. l |5.29t for a = 1 (solid line), corresponding to decoupled lay- 
ers, a = 0.5 (dashed line) and a = 0.1 (dotted line). 

z(f+ 1/2)- . = f-1/2. (5.28) 

u / > 2(1 + z Q ) 3 1 

For any value of a we obtain f{v) from Eq. ( 15.281 ), with the 
result 



2(1 -z) 



1 + z a 



where 



-I/O) 



(5.29) 



(5.30) 



The v(f) curves obtained by inverting Eq. ( 15.291 ) are shown 
in Fig. |5.8| for a few values of a. These curves resemble those 
obtained in mean-field theory and suggest the possibility of a 
velocity jump. However, before making any conclusions we 
must study the stability of these periodic orbits. An analytic 
solution of Eq. ( 15.291 ) can be obtained for a = 1/2. For a = 1, 
corresponding to rj = (decoupled particles), xq is undeter- 
mined. For any xq one recovers the single-particle result given 
in Eq. ( 15.131 ). For instance, for xq = we obtain 



/ - 1/2 
/ + 1/2 



(5.31) 
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(5.32) 



which yields Eq. ( 15.131 l. 



b. Stability. The stability of the periodic solutions found 
in the preceding paragraph can be examined from the linear 
response to a perturbation of the initial condition. Letting 
Xq + Sx, we define a Lyapunov exponent A 12 



•'•() 



x'(xq + Sx) = Xq + XSx 



(5.33) 



where x'(x ) is given by the right-hand side of Eq. (15.23b . We 
find 



A 



z 1+a {a + z a + 2/(1 + z a )) 



2fz(l + z a ) + z 1+a - az a 



(5.34) 



At the fixed point z and / are related by Eq. ( 15.291 l. Inserting 
this into Eq. (15.34b , we obtain 



A = 



z 1+a (l + z a + o(l - z)) 



z(l + z a ) - az a (l - z) 



(5.35) 



This Lyapunov exponent is plotted in Fig. l5.9l for a few values 
of a. It equals 1 for a = (corresponding to r} — > oo) and for 
a = 1 (corresponding to i) = 0). For all other values of a one 
finds A < 1 only for very small z, i.e. small v. This region cor- 
responds to the part of the v(f ) curve that has negative slope 
near the depinning threshold / = 1/2. The conclusion of this 
analysis is that this part is stable, while the portion with pos- 
itive slope is unstable. This result is somewhat surprising in 
view of the results obtained in mean-field theory. However, 
as was explained in the previous section, reversing the sign of 
the force landscape would exchange the attractive and repul- 
sive trajectories and result in the v(f) curve more similar to 
the one shown in the mean-field section. The special nature 
of the two-particle scalloped-force landscape may be related 
to the absence of zero-force saddle points which played an 
important role in the case of the smooth potential. Note that 
we have not looked for more complicated periodic solutions, 
which are difficult to rule out. 

C. Discussion of toy models 

We conclude from the previous two sections that a large vari- 
ety of behaviors can already occur with two viscously coupled 
degrees of freedom in a random-force landscape. Understand- 
ing their systematics, for instance how one evolves from the 
smooth potential to the scalloped one as more harmonics are 
included, remains to be done. In each case one must iden- 
tify the periodic trajectories and the attractor, whose structure 
may become more complex if the landscape contains more 
harmonics and more zero-force points. It is clear that an even 



12 By contrast to the Lyapunov exponent for a continuous-time flow denned 
in the previous section, this is the Lyapunov exponent for the discrete map, 
hence the transition of stable to unstable occurs as | A crosses unity. 
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FIG. 5.9: The figure shows the Lyapunov exponent A as function of 
v for a = 0.1 (dotted), a = 0.5 (dashed) and a — 1 (solid). For 
a = 1, A = 1. 



more careful systematic study is necessary when increasing 
the number of degrees of freedom within these coupled layer 
models. It is not clear at this stage whether chaotic attractors 
exist, or whether even multiple stable periodic attractors do 
coexist. 



VI. DISCUSSION 

We have studied in this paper a model of two (single- 
component N = 1) elastic layers driven over a random sub- 
strate and only coupled by a viscous coupling 7 12 , going be- 
yond mean-field theory. We have extended the functional RG 
approach which allows to describe the elastic depinning in 
each layer in the absence of an elastic coupling to the case of a 
non-zero 712. We have found that the FRG fixed point which 
describes elastic depinning is unstable to an arbitrarily weak 
viscous coupling beyond a plastic scale L p \ which diverges 
with a universal exponent as 712 — > 0. To describe the plas- 
tic physics beyond that scale we have studied the FRG to one 
loop in the moving state at non-zero velocity. We found that 
the high- velocity branch of the «(/) versus / curve terminates 
at a point v = v c where the slope is infinite. This point corre- 
sponds to a force smaller than the elastic depinning threshold, 
hence there is a range of values of / where a pinned state co- 
exists with a moving state. This dynamical hysteresis is very 
similar to the one found in mean-field theory. One could then 
conclude that the 1-loop FRG result nicely confirms the main 
features of the mean-field theory and, in addition, allows to es- 
tablish precise universal results and identify the proper length 
scales. 

This conclusion may, however, be too hurried as, surpris- 
ingly, our two-loop calculation shows some possible problems 
with this picture. The calculation is based on certain assump- 
tions (discussed in Section [IV] i.e. the neglect of violations 
of Middleton's theorem, the neglect of higher orders in the re- 
gion of small rj + ) and more work is clearly needed to ascertain 
its validity. However, as a preliminary step it indicates a new 
universality class in the case of non-periodic objects and in 
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the periodic case, a possible breakdown of the 1-loop picture 
depending on how irrelevant operators are taken into account. 
This possible alternative picture is the absence of a dynami- 
cal hysteresis and a nearly vertical «(/) curve near the elastic 
threshold. 

In an attempt to understand which effects could be missed 
by the mean field and 1-loop approaches, we have solved sim- 
ple toy models in d = 0. We found indeed that dynamical 
hysteresis may or may not be present, depending on the re- 
alization of disorder. Although in all cases considered one 
finds a periodic trajectory with non-zero velocity which sur- 
vives below the elastic (uncoupled) threshold, this trajectory 
may be attractive or repulsive depending on the disorder real- 
ization. It remains to be studied in detail how these properties 
carry to a larger number of degrees of freedom. In any case 
it cannot be assumed that a single attractive periodic attractor 
exists and a detailed study of such attractors as the number of 
particles increases must be done with care before any conclu- 
sion can be drawn. 

The particle models also show the importance of the zero- 
force points in phase space. These are couples of configura- 
tions in the two layers (uj., u 2 ) where all forces vanish. They 
are defined in the statics, hence are independent of the vis- 
cous coupling. However it is important to know their dynam- 
ical stability in presence of a viscous coupling. It is proved 
in Appendix lAl that metastable states, i.e. stable states where 
the energy has a local minimum, remain stable, i.e. dynam- 
ically attractive, at non-zero viscous coupling. This means 
that each pinned state, i.e. blocking configuration for each un- 
coupled layer, keeps a non-zero basin of attraction when the 
viscous coupling is increased from zero. In the absence of 
a viscous coupling, when upon increasing / such a pinned 
state becomes unstable, the next configuration is determined 
by the no-passing rule and Middleton's theorem, as the min- 
imum over u of all metastable configurations in the direction 
of the force. In the presence of a viscous coupling however, 
there is no guarantee that the system will not flow from there 
to a periodic orbit, resulting in a jump in the v(f) curve. Thus 
even if the metastable zero-force couples remain attractive, 
they may not be easily dynamically accessible, i.e. their basin 
of attraction may shrink and be nearly invisible in a procedure 
such as force ramping. These effects, as well as the com- 
petition between the zero-force fixed points and the periodic 
orbits, clearly remain to be studied systematically. 

It is thus a remaining challenge to understand how the FRG 
can describe the structure of such a complicated phase space 
with periodic orbits coexisting with pinned fixed points. It 
is probable that the FRG calculation presented here retains 
only averaged effects and does not adress these issues with 
sufficient accuracy. One possible geometry to study this in a 
controlled manner starting with a particle and then extending 
to manifolds is to use the drive by an harmonic well. Since 
the coupling between the layers exists only when the system 
is moving, one needs to go beyond the present calculation and 
study for instance how the avalanches in the two layers are 
correlated. 

In a broader context, one needs to justify why effort should 
be devoted to clarify the behaviour of this simple two-layer 



model with viscous coupling only. This model indeed neglects 
the competition between plastic inter-layer and elastic intra- 
layer couplings, and one should question its range of valid- 
ity. Since any small realistic interaction between two identical 
elastic layers, such as a crystal or a CDW in the absence of in- 
layer topological defects (which we have not included here) 
generates some small inter-layer commensurate coupling, one 
should ask whether our model is stable to that. This can be 
analyzed by adding a force —g\ sin (2-Kp(u], t — u^ t )), with 
p = 1, to the equation of motion ( 11.51 ) of layer 1, and its oppo- 
site for layer 2. Such a coupling generates an elastic coupling 
at small scale between the layers, and if it is relevant in the RG 
sense, at large scales as well. In the latter situation the elas- 
tic coupling sould dominate the viscous one and one expects 
that the system is described by elastic depinning. It was found 
in mean-field models that such a coupling is always relevant 
in the moving phase If such a result were general, the 
model studied here would be somewhat artificial, or describe 
only a limited range of length scales. 

It is easy to compute, to 1-loop accuracy and for any p, the 
linear eigenvalue of an infinitesimal perturbation g\ > at 
the unperturbed quasi-static depinning fixed point at v = + 
studied here. A first (and naive) calculation in the spirit of the 
statics yields: 



di gx = 2 - 4ttVA(0) 



.91 



(6.1) 



where gi is the dimensionless coupling. This is essentially the 
result obtained for the problem of a single layer in pr esence 
of both disorder and a commensurate potential 1 105] up to 
a factor of two which accounts f or the fact that the disorder 
exists in both layers, compared to Il05ll . Inserting the value of 
A(0) = e/36 at the 1-loop depinning fixed point implies that 
this coupling is always relevant near d = 4 which would seem 
to confirm the mean-field conclusion. It also yields a critical 
dimension near d = 2, i.e. 4 — d c ~ — 18/tt 2 below which the 
coupling should become irrelevant. This conclusion is correct 
for the statics, but incorrect near the depinning threshold. 

At depinning at least two new effects should be taken into 
account. First, one sees that the coupling g\ generates in 
perturbation theory a correction to the critical force, which 
amounts to adding the same force —g^ cos (2ir p{u\. t — u^ t )), 
with p = 1, to the equation of motion ( 11.51 ) of each layer. 
These two terms feed into each others and the correct lin- 
earized RG equation takes the form, to one loop: 



dfgi 
dm 



2-4ttVA(0) -27rpA'(0+) \ ( h 
27rpA'(0+) 2-4ttVA(0) J \g 2 



(6.2) 

Second, and most importantly, from the two-loop solution 
of the standard depinning fixed point [18], we know that 
A(0) does not flow to a fixed point, it always increases as 
A(0) ~ e d ~ L e . Physically, a static random force is gener- 
ated by the quenched disorder in the limit v = + . This is due 
to terms in the two-loop beta-function which account for the 
irre versib ility of depinning, and is at variance with the stat- 
ics ill 1 Oil - At depinning however, a small coupling between 
layers is always irrelevant for d < 4. This can be seen from 
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( 16.2b . since A(0) grows while A'(0 + ) converges to a O(e) 
fixed point. It justifies a posteriori the model studied in this 
paper. Of course at larger bare couplings it is likely that a 
coupled phase will arise and it would be interesting to study 
that transition ! 11 ill . 

Let us finish by recalling that one issue in the theory of plas- 
tic flow is whether one can use P(v), the distribution of time- 
averaged individual particle velocities, as a meaningful order 
parameter in the thermodynamics limit. One could then dis- 
tinguish two classes of plastic flow (i) flows with non-trivial 
P(v) (e.g. pinned particles coexisting with flowing rivers) (ii) 
flows with peaked P(v) (a delta function in the large-size 
limit) but which cannot be described by a fully elastic the- 
ory. The layered model studied here is a tractable example of 
class (ii) and requires, to exhibit a non-zero depinning thresh- 
old, elastic interactions inside the layers. Models for class 
(i) have been studied, where particles interact only through a 
hard core interaction 18211 . It would be quite interesting to find 
a tractable model which encompasses both classes and then- 
possible transitions. 
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APPENDIX A: STABILITY OF ZERO-FORCE 
FIXED POINTS 

Let us call u l x a static configuration where the force is zero, i.e. 
F z (u l (x),x) = 0, i = 1,2. The equation of motion linearized 
around the FP is: 



= M 



M = AB = - 



ill 



A -i _ I 7n 7i2 
7i2 7n 



B = 



7n #i -712-H2 

-712-ffl 7lli?2 

ffi 
-H 2 



(A.l) 
(A.2) 
(A3) 



and we are interested in the Lyapunov exponents, or relax- 
ation rates around the zero-force fixed point, i.e the eigen- 
values of the matrix M. We have introduced the Hessian 
{Hi) xx > = -V 2 x 5 xxl + V('{u l x ,x)5 XX ' in each layer, which are 
hermitian matrices. They have eigenvalues [ii. ai and eigen- 
vectors cf>i t0!i (x). In the absence of a coupling between the 
layers (712 = 0) the eigenvalues m jai are proportional to the 
Lyapunov exponents, i.e. A, iQ!j = — 7 : f 1 Vi, ai . A question is 
how they vary as the viscous coupling is increased. Note that 
one can decompose: 



(t>i, ai {x)u ai {t) 



(A.4) 



J 



and in that basis the equation of motion reads: 



E 



Tll^ai.oJ 712 f x <l>*l, ai (x)<p2, a ' 2 

r ri2l x i > 2,o, 3 ( x m,a' 1 (x) 7ll<5 



Ml,QiWai(*) 



(A.5) 



since the velocity coupling between layers is local in space, it 
becomes non-local in the eigenstates of the two Hessians. 

The matrix M has several interesting properties. Although 
it is not Hermitian, since A and B do not commute, its eigen- 
values are real. Indeed consider an eigenstate v such that 
M ■ v = Xv. This implies B ■ v = XA^ 1 ■ v, hence: 

w f • B ■ v = X v f ■ A" 1 ■ v (A.6) 

Since B and A are Hermitian (and also real symmetric) ma- 
trices, t>t • B ■ v and ■ A^ 1 ■ v are real, hence A is real. 

Consider now a bare model such that A^ 1 is strictly posi- 
tive definite with eigenvalues 7+ > 0, 7_ > 0, i.e. 7^ < r yf 1 . 
Then ( IA.6I ) implies that the sign of A is the same as the sign 
of v* ■ B ■ v. Let us consider a stable (i.e. attractive) zero- 
force point with all Hi t0li > 0, hence B is strictly negative 
definite. In its neighborhood in phase space, in the absence of 



viscous coupling between layers, the system is pinned. Since 
V' ■ B ■ v < for any non-zero v, the above property implies 
that the zero-force fixed point remains stable, i.e. all Lyapunov 
exponents remain strictly negative, as the viscous coupling be- 
tween layers 7^ < 7^ is increased, and ( |A.6t implies the 
bounds 

Mmm < _ A < AWx (A7) 

7- 7+ 

for a model with 712 < 0, and where /U m i n and /^ max are the 
smallest and largest eigen- values of —B. 

These eigenvalues of stability can be obtained exactly in the 
case where H\ and H2 commute. Then one can choose the 
same basis in both layers <pi yCC (x) = (f>2. a (x)- The Lyapunov 
exponents, i.e. the eigenvalues A in ii = Xu of (IA.3b can then 
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be organized in pairs with: 

, 7ll(Ml,a +M2,a) 



7i 2 2 ) 



± 



2(7i 2 i 

2(7i 2 !-7i 2 2 ) 



2. 2-loop order: List of diagrams 

There are seven contributions, drawn on figure uTTI Their con- 
tribution to 7 is symbolically 



(A.8) 



x4x2[a + b + c + d + c + f+g] 



(B.3) 



and one checks that as long as 7+ =711 — 712 > a stable 
FP remains stable as -f 2 2 i s increased (this holds for the two- 
particle model considered above). In general one does not 
expect Hi and H 2 to commute, since the disorders in the two 
layers are uncorrelated. For small interlayer coupling one can 
apply second-order perturbation theory: 



Ai jai 



711 



2 2 
7n - 7i2 



Ml.Ql 



(A.9) 



E 



7l2 



7ll ~ Ml.ai 



^2,02 



■|(l,ai|2,a 2 )| 2 +0(7f2) 



which always makes the smallest eigenvalue (assumed all pos- 
itive) get closer to zero, but even in the most dangerous case 
when this eigenvalue is near marginal, i.e. /Ui, ai > near 
zero, the second-order correction vanishes as /ii )Q:i — ► 0. 
Hence there is no mecanism for it to cross zero. This is not 
too surprising since the determinant of the matrix in dA.3b 



The combinatorial factor is 1/8 from the interaction, 4 from 
the time-ordering of the vertices, and an additional factor of 2 
for the symmetry of diagrams a, b, e, f and g. 

The diagrams are calculated as in 1 1811 . When expanding the 
argument of A(u xt — u xt ,), it is important to keep the index of 
the field. Only diagrams with one disorder on one layer, and 
one disorder on the other layer can give rise to a contribution 
to 712, which will be the new feature found below. 

3. Expressions for the diagrams 

The first combination is 

a + g= -A"{0+) 2 lf (B.4) 



as before, since the free integration kills the inter-layer term. 
In the following we give corrections proportional to u\. The 
index i runs over both layers. Integrations over momenta and 
time are not written. 

nt of the matrix in dA.3b 

cannot change sign as 7l 2 2 < 7^ is increased. Since the b + c + d = Ui ^Ru{qi, h)R n (q2, h)Ru{qi, £3) 



eigenvalues remain real (as shown above) they cannot con- 
tinuously change sign. Hence, as above, we conclude that a 
stable zero-force fixed-point remains stable. Both the size (in 
phase space) of its basin of attraction (pinned phase) and the 
Lyapunov exponent may decrease with increasing interlayer 
coupling, but they do not cross zero. 

Finally note that if e.g. Hi has a marginal direction, i.e. 
Hi-v 1 = 0, then v = (v 1 , 0) is an eigenvector of M with zero 
eigenvalue. Hence a marginal direction remains marginal. 



APPENDIX B: 2-LOOP CALCULATIONS 

In this appendix we derive the FRG equations up to two loops 
using the method of ref. ITU]. All calculations are done at 
zero velocity, at the depinning transition. All static quantities 
like the disorder correlator are the same as for the standard 
depinning transition, and we refer to lfl8ll for details. Here 
we only calculate the corrections to friction, i.e. corrections to 
711 and 712. 

1. 1 -loop order 

There are no corrections to 712 at 1-loop order, since there 
exists only a single vertex, thus one cannot get a term of the 
form J u 2 v}. For this, one needs (at least) 2 loops. Therefore: 



x [|i3-ti|-|i3 + *2-ti]A"'(0)A'(0)^ 

(B.5) 

= -^2Rii{qi,ti)Rii{q2,t 2 )Ri l {qs,h)\U - t 2 \ 

i 

xA'^OjA'tOjii 1 (B.6) 



f : 



-2A"'(0 + )A'(0 + )/a - 2A"(0+) 2 /a (B.7) 



Integrating over times yields the diagrams presented in the 
next two subsections. They involve the following non-trivial 



11 * 





Sjl\ oop = -A"(0)/i7ii 



(B.l) 






h 



p (p 2 + m 2 ) 2 



(B.2) 



FIG. B.l: 2-loop dynamical diagrams correcting the friction. They 
all have multiplicity 8 except (c) and (d) which have multiplicity 4. 
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momentum integrals: 



L, := 



Ia := 



(qf + m 2 )(q 2 + m 2 ) 2 (g 2 + q 2 + 2m 2 ) 



1 l-21n2 



2e 2 4e 
d d qi d d q 2 1 



(e/i) 2 + finite 



(B.8) 



1 



{2w) d (2n) d q\ + m 2 g 2 + m 2 ((gi+9 2 ) 2 + m 2 ) 2 

(B.9) 



They are calculated in lfl8ll . 



1 



9192 («? + m 2 )(q 2 + m 2 f{ql + q\ + 2m 2 ) 



as can be seen by symmetrizing in qi and q 2 . 

4. Corrections to 712 



hi 



a+g 







h 2 

2 h ' 
(B.10) 



(B.ll) 



h h it 



1I2 ((ll + 4) ill + fe 2 - ql) 1I2) A'(0+)A"'(0+) 



9192 2q 2 q 2 (q 2 + q 2 )j 11 {{q 2 



llflli 



{q\ - 9|) 2 7i 2 2 ) 



r 



(B.12) 



In principle, ( IB. 12b should be written with massive propaga- 
tors. We have put m — for notational compactness. It is 
easy to see that for generic values of 711 and 712 (1B.121 > has 
no subdivergence for either qi — > 00 or q 2 —►00, only if they 
become large together. Properly regularized, it therefore has 
only a single pole in e, and this pole is universal, i.e. inde- 
pendent of the regularization scheme. For such an integral, 
which moreover is homogenous in qi and q 2 , the pole can be 
expressed as 



f(ql,q 2 ,m 2 ) 



(eh? 



d(q 2 2 )f(l,q 2 ,0)q 2 +O(e ) 



(B.13) 

This is proven using conformal m appings of the d ifferent sec- 
tors, and was established in Col[!o3,[!ol[l09]. Accepting 
that the integral is indeed universal, a quick way of deriving 
( IB. 13b is as follows 



f(ql,q 2 ,m 2 ) 



f{qt,q z 2l 0)~F- 

lqi<AJq 2 e 

(B.14) 

where the "w" indicates up to terms of O(e). To obtain the 
residue F, we derive w.r.t. A, and then set A = 1: 



F^Sr 



/(<h 2 = 1,^,0) 



(B.15) 



Since the integral is finite, we can take the limit of e — > or 
d — ► 4. This gives the result ( IB. 13b up to an overall normal- 
ization, which is also easily checked. 

It is now straightforward to integrate ( IB. 121 ) using ( IB. 131 ): 



h 



7l2 A'(0+)A( 3 )(0) 



12 



4e 



In 



7n + 712 



711 — 7i2 



(B.16) 



Note that physically one has to restrict to 7+ > and 7_ > 
0, thus all results are to be taken in this domain only. The 
absolute-value therefore represents nothing but a notational 
commodity. 

The other two diagrams are trivial: 



Sj c 12 = 

Hi = • 

5. Corrections to 711 



(B.17) 
(B.18) 



Grouping diagrams, which partially cancel, we find for the 
corrections to 711: 



«57 1 T g =7iiA"(0+) 2 / 1 2 



(B.19) 



c b+c+d 
<*7ll 



(2 (ql + g 2 2 ) 2 7ii + (-H + 3g 2 2 9i 2 + 4) ihlli + 1 2 2 (ql - ql) 7i 4 2 ) A'(0+)A'"(0+) 

1^ ^qlqi (ql + ql) m ((ql + q\f iii - (ql ~ ql? 7 2 2 ) 

712=0 r 7ll A'(0+)A"'(0+) 



274T7^-2^ (B.20) 



9192 



qm (ql + ql) 
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^7ii~ c+d i s a little bit more complicated to calculate analytically, since it has a subdivergence, which has to be subtracted if 
we want to use the magic relation ( IB. 13b . We observe that subtracting the term at 712 = 0, the diagram has no longer a 
subdivergence. In order to proceed, one then uses ( IB.13l l. and performs a partial fraction decomposition (the variable is g|) of 
the remaining term, leading to integrals known by Mathematica. The final result is 



b+c+d _ 7ll A'(0+)A"'(0+) 2 7l2 A'(0+)A"'(0+) 



7n + 7i2 



7n — 7i2 



(B.21) 



The next diagram is 

*7i e i = 



711 


{(4 + 4) 2 (4 + 4) 1I1 - (4 - 244 - H4 H4 + 4) ih) 


A'(0+)A"'(0+) 


H44 (4 + 4) ((4 + 4? 7?i - (4 - 4) 2 7? 2 ) 





(B.22) 



We have used the abbreviations 173 := qi + q 2 , and q 3 := | ^3 1 . 
Again, this diagram has a subdivergence (double pole), which 
we want to subtract. Let us again try the term at j 12 — > 0: 



Slli 



7w-o f ( (? | + (? 3 4 )7 11 A'(0+)A'"(0+) 



2ql44(4 + 4) 
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,2„4 („2 



< qiq2 44(4 + 4) 

S^fi can be rewritten as: 



= A'(0)A"'(0) 



9192 



(4 + 4) in 

H44(4 + 4) 
27n7 2 2 



4 (4 + 4f 1I1 - 4 [4 - 4? (4 + 4) ih 



A'(0)A"'(0) 



9i 92 



711 



44(4 + 4) 



j 



27ii7i 2 2 



4 (4 + 4) ((4 + 4? ill (4 4 f ih) 

(B.24) 



(B 23) Tn e ^ ast i nte g ra l i s 



» 4 (4 + 4) ((4 + 4? ill - (4 - 4? ih) 



log 



711+712 
7ll — 712 



log 



7l2 

7?j~ 



4e7n7i2 



4e7i 2 2 



(B.25) 



One way to prove this is as follows. Introduce Schwinger-parameters to write the l.h.s. of (IB.25b as 

/ / e -Si(92+93) 2 -(92+93) s 2-S3(7ll92+7l292+937ll-937l2)-S4(7ll92-7l292+937ll+937l2) e -S 2 (B 26) 

Jqi,q 2 J si>0 : s 2 >0,s a >0,si>0 

where we have introduced a mass for s 2 only (using again universality of the leading pole in e). Then integrate over the q^'s: 

-d/2 



I 



si>0,s 2 >0,S3>0,s 4 >0 



(- (S3 ~ S4) 2 7i2 + s 2 (2si + s 2 ) + (s 3 + s 4 ) 711 (2 (si + s 2 ) + (s 3 + s 4 ) 7u ) 

I 



(B.27) 



Rescale all Sj with i ^ 2 by s 2 and integrate over s 2 . Then 
go to new variables S3 — > (s + t)/2, s 4 — ► (s — t)/2. Our 
integral becomes 



s>0 



dt 



Sl >o 2 ((s 7ll + 1) (2sj + S711 + 1) - t 2 7l 2 2 ) 

(B.28) 

The result can be simplified to (IB. 25b . A tricky point are logs 
halfway. Expanding dB.281 > in 712, we circumvent the prob- 
lem and can check the first terms of the Taylor series. The 



complete result for is (up to finite terms) 



d^n = A'(0)A"'(0) 



7iiA 



712 
" 2e 



log 



711 + 712 



7n — 7i2 
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log 
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The final diagram is 



(B.29) 
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f 7ll ((g 2 + gl) A'(0+)A"'(0+)g 2 + g 2 (g 2 2 + gf) A"(0) 2 ) y A'(0+)A"'(0+) + A"(0) 2 

711 = / T^v = 711 / 

Jqiq-2 </l ( /2 ( /3 Jqi 

= 2711 [A'(0 + )A"'(0 + ) + A"(0) 2 ] Ia ■ (B.30) 
This gives the flow equations given in the main text. 
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